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ABSTRACT 

After  a  brief  introduction  on  the  parameter  plane  methods,  the  be- 
havior of  simple  systems  with  time  delay  has  been  initially  studied  by 
the  root-loci  method  and  then  the  state  of  the  art  of  parameter  plane 
methods  with  time  delay  has  been  presented. 

A  variation  of  the  solution  in  the  state  of  the  art  presentation 
has  been  derived  and  an  appropriate  digital  computer  program  has  been 
constructed.   As  an  extension  of  studies  in  the  parameter  plane,  a  three 
variable  parameter  system  has  been  approached  by  using  basic  descriptive 
geometry  properties  and  simple  illustrative  examples  have  been  worked 
out  in  the  parameter  space. 

Finally,  a  real  life  system  with  time  delay  has  been  analysed  and 
redesigned  by  parameter  plane  methods. 
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I.   BASIC  PARAMETER  PLANE  TECHNIQUES 

A.  GENERAL 

Parameter  plane  techniques  are  algebraic  methods  designed  to  study 
the  behavior  of  linear  systems  through  their  characteristic  equation, 
when  two  of  their  parameters  are  adjustable.   So  in  fact  they  are  an  ex- 
tension of  the  Root-Locus  technique  where  only  one  of  the  parameters  of 

# 

the  system  is  adjusted  to  meet  any  desired  performance  requirements. 

B.  SILJAK'S  APPROACH 

Parameter  plane  techniques  have  been  at  first  studied  by  Mitrovic 
who  considered  that  the  parameter  variations  directly  affected  two  of  the 
coefficients  of  the  system's  characteristic  equation. 

Then  Siljak  extended  Mitrovic ' s  work  by  including  variations  on  the 

system's  parameters  "or"  and  "|3"  which  appeared  linearly  in  one  or  more 

of  the  coefficients  of  the  characteristic  equation,  which  was  of  the 

form: 

n 
F(s)  =  £  aksk  =  0  (1-1) 

k=0 

where  a,     =   bo-  +  c,  p  +  d,  (I-la) 

with  b.,  c.,  d.  constant  and  real 
1'   1'   1 

and  a, (3  variable  and  real. 
Curves  obtained  in  the  a, (3  plane  are  called  parameter  plane  curves 
and  they  provide  information  about  the  roots  of  the  system's  character- 
istic equation  in  a  similar  way  to  root  loci. 
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Siljak's  work  can  be  summarized  in  its  goals  as  follows: 

First  make  an  appropriate  substitution  for  s  in  equation  (1-1).   This 
substitution  must  be  in  terms  of  F    and  uu  ,  or  a  and  uu,  or  in  general  in 
terms  of  quantities  one  wants  to  study. 

Second  step  is  to  separate  the  result  of  equation  (1-1)  after  the 
substitution  is  made,  into  real  and  imaginary  parts  thereby  obtaining 
two  equations  in  o/,f3,  £,  uu  and  the  constants  of  the  system. 

Finally,  one  solves  for  a,  (3  in  terms  of  £  and  uu.   If,  for  example, 
one  fixes  £  =  T   and  varies  uu  over  the  range  0  <  uu  <  »,  one  will  obtain 
in  the  a-p  plane  the  necessary  relation  between  a   and  £5  in  order  to  have 
a  complex  root  pair  of  damping  ratio  £  =  £   (a  constant  £  curve).   If, 
on  the  other  hand,  one  fixes  uu  =  uu  and  varies  £  over  the  range  -1  <  £  <  1, 
one  will  obtain  a  constant  omega  curve  denoting  the  necessary  a-f3  relation 
to  provide  a  complex  root  pair  of  natural  frequency  UU  =  uu  .   Siljak  used 
the  following  substitution  in  order  to  study  constant  £  and  constant  uu 
curves : 


k    k 

s  =  uu 


_Tk  «)  +  2  J   1-C2  V<)]  (1-2) 


where  U,  and  T,  are  Chebyshev  functions  of  the  first  and 
second  kind  respectively,  given  by  the  following  recursive  relations  and 
initial  conditions: 

Tk+1"    2ftk(O   +  Tk.l(O   =   0 

(1-3) 
Uk+l"    2?W   +Uk-l<t>   =   ° 

where 

T0(C)   =   1  U   <£)  =  0 

°  and  °  (1-4) 

tx<E)  =  z  ux(r)  =  1 
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It  is  to  be  noted  that: 

\(-0  =  (-Dk  \(0 

k+1  (I"5) 

uk(-C)  =  (-D   uk(g) 

In  the  second  step  after  substituting  equation  (1-2),  [considering 

also  the  results  of  equations  (1-3)  to  (1-5)],  into  equation  (1-1),  Siljak 

obtained: 

n 

V  (-l)k  akwk  Uk-1(C)  =  0   and 


k=0 
n 


I  <■-»*  v"  uk  «)    =  o 


So  finally: 


BjCV  +  C.p  =  -D, 


B2<*  +  C2P  =  -D2 


from  which: 


(1-6) 


k=0 
where  a,  =  b.ot  +   c,  p  +  d  (I- la) 


(1-7) 


Cl  D2  -  C2  Dl  p  =  B2  Dl  -  Bl  D2 

Bl  C2  '  B2  Cl  Bl  C2  "  B2  Cl 

where: 

n  n 

bi  -  l  <-x>  V\-i  B2  - 1  (-1)  Vkuk 

k=0  k=0 

n  n 

k=0  k=0 

n  n 

"i'I^'Wi  D2 = I  <-«kv\ 

k=0  k=0 

The  study  for  conditions  in  the  o>p  plane  giving  a  real  root  at 
s  =  a,  was  much  easier.   In  fact,  by  substituting  s  =  a  in  equation  (1-1), 
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Siljak   obtained 
n 

£  akak  =   0  (1-9) 

k=0 

which  by  using  (I- la),  led  to: 

n         n        n 

V    k    V    k   V    k 
a)   b,  a     +   B)  c.  a  +  )  d,  a  =0 

L  k      Z_i  k     /_>  k 

k=0       k=0      k=0 
which  is  the  equation  of  a  straight  line  in  the  Qf-B  plane,  for  given  "a". 
Details  on  this  work  can  be  found  in  reference  [lj. 

C.  HOLLISTER'S  WORK 

In  section  4  of  reference  [2]  Hollister  obtained  parameter  plane 
equations  for  a  system  with  characteristic  equation  of  the  form  of  equa- 
tion (1-9),  with: 

ak   =  \A  +  Cl-B  +  dkAB  +  Ac  C1"10) 

His  solution  for  constant  "a;  "  and  "Z"  curves  was  obtained  by  simul- 

n 

taneously  solving  two  equations  of  the  form: 

B.A  +  C  B  +  D  AB  +  F   =  0 
i-  i-  L  L  (i-U) 

B2A  +  C2B  +  D2AB  +  F2  =  0 

where  B.,  C . ,  D.,  F.,  with  i  =  [1,2],  are  functions  of  "cu  "  and  "rn,  as 
1111  '   '  n 

well  as  of  the  first  and  second  kind  Tchebychef  functions  T,  (C)  and  U,  (C) 

D.  PRACTICAL  CONSIDERATIONS  (AN  INDICATION) 

Although  the  concepts  governing  analysis  and  synthesis  of  feedback 
control  systems,  once  established  are  quite  straightforward,   there  is 
some  difficulty  in  making  an  appropriate  selection  of  the  range  of  A  and 
B  parameter  values  to  be  used  in  the  digital  computer  program. 

Some  practical  ways  for  defining  a  set  of  approximate  ranges  of  the 
parameter  values  for  use  in  the  first  run  of  a  parameter  plane  computer 
program,  are  presented  in  section  V. 
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II.   SYSTEMS  WITH  TIME  DELAY  --  A  RESUME 

The  purpose  of  this  section  is  two- fold: 

First:   To  introduce  the  reader  to  the  State  of  the  Art  of  analysis 
for  systems  with  transport  lag. 

Second:   To  define  the  objectives  of  the  work  presented  in  this 
thesis. 

Since  stability  analysis  will  be  treated  separately  in  a  later  sec- 
tion, items  concerning  this  subject  are  not  included  here. 

A.    THE  TRANSPORT  LAG  -  AN  INTRODUCTION 

1.    The  existence  of  a  time  delay  effect  associated  with  any  sys- 
tem, linear  or  non-linear,  is  inevitable  since  the  working  signal  (com- 
mand or  information),  needs  a  certain  period  of  time  to  travel  through 
and  reach  the  output.   So  if  a  signal  r(t)  is  applied  in  the  input  of 
one  of  the  most  elementary  systems,  as  shown  in  figure  2-1,  the  output 
c(t)  will  be  in  general  equal  to: 

c(t)  =  Kr(t  -  t)  (II-l) 

where  T  is  a  time  interval  associated  with  the  equivalent  length  of  the 
system  "£"  and  the  average  speed  u  of  the  signal  through  the  device. 
Specifically: 

T  =_J_  (II-la) 

u 

There  are  of  course  cases  when  the  dimensions  of  the  system  and  the 

speed  of  the  signal  are  such  to  imply  a  negligible  time  delay  effect  as 

for  example  in  the  actual  lumped  A.C.  circuitry  where  the  wavelength  of 

the  signal  is  large  compared  with  the  dimensions  of  the  circuit-system. 
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But  there  are  also  other  cases  when  the  time  delay  effect  is  import- 
ant and  should  be  taken  into  consideration. 

In  such  a  case  one  can  define  a  system  in  the  s-domainj equivalent 
to  the  one  of  figure  2-1  by  Laplace- trans  forming  equation  (II-l): 

C(s)  =  K  (R(s)e"TS)  =  (Ke"TS)R(s) 

C(s)       -TS  ,__  _. 

°r   R(s)  =  Ke  (II'2) 

2.    After  this  presentation  one  can  define  the  time-delay  effect  as 
"the  physical  phenomenon  of  finite  time  spent  by  a  signal  travelling 
through  a  finite  system."   For  each  combination  of  signal  and  system 
there  is  assiciated  a  time  delay  constant  "t"  which  denotes  the  time 
necessary  for  this  travel. 

B.    THE  EFFECTS  OF  TRANSPORT  LAG  ON  A  SIMPLE  LINEAR  SYSTEM 

Before  the  State  of  the  Art  in  analyzing  systems  with  transport 
lag  in  the  parameter  plane  is  discussed,  it  is  interesting  to  observe 
the  results  of  this  time  delay  in  a  simple  system  by  use  of  conventional 
methods  of  analysis. 

1.  Consider,  for  example,  the  system  of  figure  2-la,  with  trans- 
fer functions  as  follows: 

a.  Open  loop:   GH  =g/s+1\ 

b.  Closed  loop:  — — -  =  — x 

'   R      S   +  S  +  K 

From  the  root  locus  analysis  of  this  system  as  shown  in  figure  2-lb,  a 

pair  of  roots  at  Z  =  .707  and  U)  =  .707  is  derived  for  K  =  \. 

n 

2.  In  figure  2-2  the  same  system  is  shown  with  the  exception  that 
a  time  delay  is  now  present  in  the  forward  path  and  the  transfer  function 
changes  to: 
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,       _u     K -ts  (II-3a) 

a.  Open  Loop:   GH  =  g/s+]\  e 

"TS 

b.  Closed  Loop:   |  =  ■  ■_  K  e (II-3b) 

R   S  +  S  +  Ke"TS 

In  trying  to  apply  the  root-locus  analysis  technique  one  is  faced  with 
the  first  difficulty  since, 

K 


GH  = 


S(S+l)eTs 


Ts 
has  an  infinite  number  of  roots,  as  the  series  expansion  of  e  '  implies. 

3.    There  are  a  number  of  approximation  techniques  one  can  apply  to 

overcome  this  difficulty  as  mentioned  below: 

a.  For  small  T  and  low  frequency, 

esT  =  1  +  sT  (II-4) 

b.  From   the    fact    that 

e   =    lim  I  1   +  —    j 
.     \  v   / 

V-f00 

V 
A  '  -«*      T        (        l 

and  e   =    lim 


1   +  cv 
v~-><*  — 

v 


one   can   approximate: 


-sT 
e 


(      1    \V 
=  ( ~  1       if  v  is  sufficiently  large.   In 


1  +~ 

v 

order  to  have  an  intuitive  feeling  of  how  large  v  should  be  for  accept- 
able results,  note  that  (1  +  1/10)    =  2.6  producing  a  deviation  of  4.4% 
from  e  =  2.72.   So,  any  v  <  4  will  produce  too  crude  an  approximation 
since  (1  +  \)k  =   2.44,  and  2.44  r  2.72  =  9%. 

c.    In  the  case  where  both  the  frequency  and  T  are  too  large 

-Ts 
to  neglect  in  the  expansion  of  e    ,  it  is  possible  to  approximate  the 
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time  delay  term  by  the  ratio  of  two  polynomials  in  "s"  yielding  to  a  cer- 
tain order  of  terms  identical  with  those  contained  in  the  series.   That 


is, 

N  U(+Ts) 

e'Ts  =  F    (T.s)  =  -^- 

e       a,bU  S;    D    (+Ts) 

3  *  D 

where  "a"  and  "b"  are  the  order  of  the  numerator  and  denominator  poly- 
nomials.  So  since 

2      3       4 
-Ts    .    _  ^T  2   T   3   T   4 

e    =  1  -  Ts  +  -^  •  s  -  -7   *  s   +0-r  •  s   -... 
2      6      24 

if  third  order  terms  can  be  neglected,  it  can  be  shown  by  straightforward 
division  and  a  series  of  comparisons  with  the  series  terms  that  a  conven- 
ient approximation  is  the  following: 

-Ts        1+^s+ApS2  1    -    3/4(Ts)+  T2/4(s2) 

6  1  +m>is  L  +  T/4(s) 

These  are  the  basic  features  of  the  so-called  "Pade  Approximation  Tech- 
nique."  In  appendix  A  the  most  often  used  approximations  are  listed 
(degree  of  higher  order  polynomial  equal  to  12).   In  this  introductory 
study  of  time  delay  effects  the  Pade  approximation  will  be  selected  as 
the  most  accurate.   From  appendix  A  then, 

-(Ts)             N4  4(Ts) 
e  K      }   =   F    (Ts)  =  — -^ = 

4>4^      D4,4(Ts) 

2  2      3  3    4  4 
1680  -  840Ts  +  180T  s   -  20T  s   +  T s 


1680  +  840Ts  +  180T2s2  +  20T3s3  +  T4S4 
and  by  use  of  this  approximation  the  closed  loop  transfer  function  of 
the  system  becomes: 

C        K#N4  4<Ts) 

C     fbit (H-5) 

(s  +s)04  4(Ts)  +  K-N4  4(Ts) 
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4.  The  root  loci  of  this  system  have  been  plotted  for  four  different 
values  of  T  by  use  of  a  digital  computer,  as  the  following  table  indicates: 

TABLE  2-1 

T  =  0  Figure  2- lb 

T  =  .1  Figure  2-2a,  2-3a 

T  =  .2  Figure  2-2b,  2-3b 

T  =  .3  Figure  2-2c,  2-3c 

T  =1.0  Figure  2-2d,  2-3d 

It  should  be  noted  that  the  Pade  approximation  introduces  in  general, 
roots  in  the  system  which  all-together  do  not  disturb  seriously  the  dom- 
inant roots  near  the  origin.   Data  on  the  characteristic  equation  result- 
ing for  each  of  the  above  "T's"  are  shown  in  Table  A2  of  Appendix  A. 

5 .  Partial  Conclusions 

a.   A  careful  observation  of  figures  2-lb,  2-2,  2-3,  and  Table 
2-1  shows  some  of  the  effects  of  the  existence  of  transport  lag  in  the 
linear  system  under  consideration,  namely: 

(1)  With  the  insertion  of  time  delay  in  the  system  the 
segment  of  the  root  locus  parallel  to  the  jjj  axis  doesn't  exist,  but  with 
increasing  time  delay  constant  the  corresponding  segment  is  bent  analog- 
ously toward  the  juu  axis.   In  other  words  the  system  becomes  unstable. 

(2)  As  a  consequence  the  K    of  the  system  which  main- 

n  max 

tains  absolute  stability  is  decreased  with  increasing  time  delay. 

(3)  It  is  also  obvious  (Table  2-Z) ,  that  the  relative 
stability  is  reduced  with  increasing  time  delay  constant. 
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6.    General  Conclusions 

From  this  discussion,  the  main  facts  concerning  the  insertion  of 
a  time  delay  in  a  linear  system  are  apparent: 

a.  Any  system  carries  a  time-delay  effect  with  it.   Therefore, 
no  system  can  be  thought  of  as  absolutely  stable.   Nevertheless,  systems 
with  very  small  time  delay  deviate  so  little  from  the  ideal  case,  that 
they  can  be  thought  of  as  delay- free. 

b.  By  considering  Figures  2-3  we  see  that  when  approximating 
the  time  delay  term  by  the  Pade   technique,  poles  and  zeroes  are  intro- 
duced at  finite  UU  ,  in  the  "s"  domain.   This  is  true  for  any  other  approx- 
imation. 

c.  For  T  increasing,  the  relative  stability  of  the  system 
decreases  and  the  root  locus  suffers  an  increasing  deformation  towards 
the  +o~  direction. 
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C.    PARAMETER  PLANE  METHODS  FOR  SYSTEMS  WITH  TIME  DELAY  (STATE  OF  THE  ART) 
1.    The  closed  loop  gain  of  the  system  shown  in  figure  2-4a  is 


R 


(s)   = 


N(s) 


D(s)eTs  +  N(s) 


with 


n 

\  k 


V 


(H-6) 


D(s)   'i   )    d   s"        and        N(s)   -   ) 
k=0  k=0 

So   the  characteristic   equation   is 


n  J  n 

Ts  V         k       V         k       V  k        _ 

e     L\s    +L\S    =  L\(s)s    =  ° 

k=0  k=0  k=0 

Ts 
where  a,  (s)  is  a  function  of  e   ,  and  contains  the  variable  parameters  A,B 

of  the  system.   So  if  only  linear  combinations  of  a  and  b  are  considered: 

a^s)  =  A'ok  +  acke  s  +  BdR  +  be^  S  +  fR  +  gRe  S     (H-7) 

2.    Intending  to  derive  constant  'T"  and  MUU  "  equations  for  the 

^        n 

system,  one  must  substitute  s  in  the  characteristic  equation  in  terms  of 
£  and  uu  ,  that  is  by: 


s  =  -Cx  +  j'ju  V  1  -  r 


(See  figure  2-5). 


Therefore,  it  follows  (Siljak)  that 


k     k 
s   =0) 


T  K)  +  jV  i  -  XT  u. .  (-C) 


w 


ith  T,  (-£)    and   U,  (-£)   being   Chebyshev   functions   of   the    first   and    second 

kind   defined  by 

\(V   =   cos(Kcos_1C)    and   U   (£)    =   sln(Kco^  C) 

sin(cos        £) 

(see  section  I-B  for  the  recursion  relations),  and  note  that 
\(-0  -  (-1)\(C)  and  Uk(-C)  =  (-l)k+1Uk(C). 
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3.    After  the  substitution  for  s  in  equation  (II-6)  is  made,  and 
by  using  the  information  in  the  previous  paragraph  and  equation  (II-7), 
the  following  result  is  obtained  after  separating  reals  from  imaginaries: 

[AB1(uun,C)  +  BC^lD^C)  +  D1(«Jn»C>]  + 
j  [AB2(uun,C)  +  BC2(Xn,D  +  D2('JJn^)]  =  °  (II~8) 

where 


Bl   =I'JJnk[("1)kTk(C)bk-H:"0Ck(("1)kTk(C)cOSe    ~(-Vk+lJl<2   uk(C)sin9 


k=0 
n 


Cl  =IU3nk[(-1>\^)dk-^"0ek((-l)\(C)cos9    -(-Dk+7l-C2   Uk(C)sine)_ 


k=0 
n 


and 


D.    =   >   tt) 

1       Li 

k=0 


n 


^mnk  [(-l)kTk(C)fk-f€-0gk((-l)kTk(C)cos9    -(-l)k+iyi-C2   Uk(C)sine)] 


\;«  kr^M 


■2.    ,    "0, 


.k+L 


/,    Jl 


k. 


B2  =A    l(4)    X^V^V     ckV(_1)     \<ZNl<  cosO+C-n^CC)  sinQ;_ 

k=0 
n 


2       L>  n 
k=0 


J-l)k+1Uk(C)^l-C2dk+e"0ek((-l)k+1Uk(C)/l-C2cose  +  (-l)kTk(C)    sine)] 


D2   "  1\ 
k=0 


(-Uk+V(CxA<2fu+€"0gJ(-l)k+V  (CxA-C2cos9  +  (-l)kT.  (C)    sin9 


|Note   the   identity     esT  =  e_0(cos9   +  jsin9)    for  0  =  £'u   T  and  9    =  uu  Vl-^      T| 

4.  The  solution   for    A  and    B  easily   follows    to  be: 

C.D?    -    DC                      D  B?    -    D„B, 
A  =  s ;      B  =  S ;      A  -   BlC2    -   ClB2 

5.  This  derivation  is  the  one  followed  by  L.  Eisenberg  [4].   A 
derivation  of  parameter  plane  equations  for  a  system  with  time  delay  and 
two  variable  parameters  included  in  non-linear  combinations,  will  follow. 
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6.    System  with  the  Product  AB  Present  in  its  Characteristic  Equation. 
Consider  the  system  in  figure  2-4b  with  closed  loop  gain 

C(s)  .  G[s) 

RK    J  Ts 

e   +  G(s)H(s) 

For 

G  (S)  H^(s) 

G(s)  =7rrT   and   H(s)  =  — 777 
GD(s)  ^^^ 

the  characteristic  equation  will  be: 

N 
F(s)    =   GD(s)HD(s)eTs   +  GN(s)HN(s)    =  £  a^  =   0  (II-9) 

k=0 


where 


ak =  A(\ +  vTS)  +  B(dk +  vTS)  +  AB(hk +  vTS)  +  (fk  +  gkeTS) 

Ts 
An  expression  for  e    can  be  derived  by  using 


/   2 
s  n    J  n  V 


and  calling 


So 


;nT  =  0  and  ^  A 


Ts    -0   j8     -0,  '     .   .  on 

e    =e^eJ   =  e  ^(cos9  +  j  sin9)  . 

k 
The  expression  for  s   as  a  function  of  <D  ,  £,  and  the  Chebyshev 

functions  are  also  available.   So  the  characteristic  equation  can  be  ex- 
pressed as  a  function  of  f    and  U)   as  follows: 

*  n 

N 


F«'V  ■  ZA<-Vk  [\«>  -  V1"?2  °k«)]  ("-10) 


where 


k=0 


a^  =  A(b     +  c  e'^cosS)   +  B(d      +   e  e~0cos9)   +    (f  +g  e"^cos9) 
+  AB(h-f   i  e"^cos8)   +  je  ~^sin9  (Ac  +  Be  +  g  +ABi    )  =    aR+Ja! 
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Next,  by  defining  £3  =  T.  (£)  and  R  =  Jl-C   Ui,(C)>  tne  characteristic 
equation  can  be  expressed  as: 

n 

F<C.»n)  =  £<-«>/  [vr  +  Vi  +  J(aiWi>]  ■  ° 

k=0 
which  after  separating  into  reals  and  imaginaries  leads  to: 


and 


Re 


Im 


[f(C»V_ 


N 

k=0 

N 

k=0 


And    finally,    by   expanding   terms    inside  parenthesis,    this   set   of  equations 
can  take   the    form: 

AB,    +   BCL  +  Dl  +   ABEj^   =   0 


AB2   +  BC2   +  D2   +  ABE2   =   0 


where: 


B 


i  "h-"A\\<-0+  v-*(Tk(?)cose  +  uk«:)/i-c2  sine)] 


k=0 

N 

1        L-i        n 
k=0 

N 


dkTk(C)+  eke"0(Tk(C)cos9    +  Uk(Cy'l-C2   sine)_ 


Dl   =l^AV^)+  Zke-\\(OcosB    +  Uk(C)/l-C2   sinejj 


k=0 

N 


El   =I(^n)k[hkTk^)+   V'0(Tk^)cOs6    +  "k^V1-^    sine); 


k=0 

N 


B2  -l<'"Ji-\\Wy/l'Z2  +   V-0(Tk(C)sine    -    ^jl-^   cose)_ 


k=0 
N 


C2=l^nA-\\^M2  +  ^'\\(OsinB    -    \QL\/l<2   cose)] 


k=0 
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N 


»i  " If-Vl" Ww/l-C2  +  ske"0(Tk«)Sme  -  \(oJi<2  cose)] 


k=0 

N 


E2  -  l^nA-\\^M2  +   V"*(lk<Osine  -  Uk(o7l-C2  cose)] 
k=0 

The  selection  then  of  (11-11)  is  a  matter  of  algebraic  manipulation. 
7.    The  derivation  in  paragraph  6  above  is  the  one  followed  by 
S.E.  Lamberski  [5], 

C.    THE  OBJECTIVES  OF  THIS  WORK 

1.  What  has  been  discussed  up  to  this  section  was  a  review  for 
people  familiar  with  control  theory  and  hopefully  an  introduction  to  the 
state  of  the  art  to  others „ 

2.  The  sections  to  follow  contain  material  associated  with  the 
fulfillment  of  objectives  as  follows: 

a.  Theoretical  study  of   general  systems  with  three  variable 
parameters  and  with  time  delay.   This  is  shown  in  Section  V. 

b.  Study  of  actual  engineering  problems  involving  two  or 
three  variable  parameters  in  systems  with  time  delay. 

c.  (1)   Since  for  such  an  engineering  study  to  be  complete, 
a  stability  analysis  is  needed  and  hence  Section  IV  is  devoted  to  pre- 
senting the  state  of  the  art. 

(2)   In  Section  III  parameter  plane  equations  are  derived 
for  a  general  system  with  two  variable  parameters  and  time  delay.   The 
reason  that  this  derivation  has  been  included  is  that  it  shows  an  inter- 
esting way  to  deal  with  the  time  delay  difficulty.   In  all  of  its  other 
aspects  it  is  the  same  as  the  methods  mentioned  in  this  section,  since 
they  involve  the  already  familiar  algebraic  manipulation  of  Siljak. 
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III.   SYSTEMS  WITH  TWO  VARIABLE  PARAMETERS  AND  TIME  DELAY 

The  purpose  of  this  section  is  to  derive  parameter  plane  equations 
for  a  general  system  with  time  delay  and  two  variable  parameters.   The 
basic  parameter  plane  method  is  used,  but  it  is  different  in  the  manner 
the  time  delay  effect  is  treated.   This  approach  will  require  a  shorter 
computer  program  to  implement. 

Since 'the   investigation  of  this  kind  of  a  system  started  with  the 
analysis  of  a  specific  actual  problem,  equations  for  that  kind  of  a  pro- 
blem are  first  derived.   Extension  to  the  general  case  is  then  covered 
in  two  further  steps,  by: 

first,  considering  other  possible  real  life  system  configurations 
and  systematically  filling  in  possible  terms  for  a  more  general 
representation  of  all  of  them  by  one  equation,  which  is  true 
under  certain  conditions, 
and       second,  removing  the  imposed  conditions  and  thus  getting  the 

general  parameter  plane  equation  for  a  system  with  two  variable 
parameters  and  time  delay. 

A.    A  SYSTEM  WITH  TIME  DELAY  IN  THE  FORWARD  PATH,  AND  THE  VARIABLE 
PARAMETERS  IN  THE  FEEDBACK. 

In  this  section,  parameter  plane  equations  for  a  system  with  time 
delay  in  the  forward  path  and  the  variable  parameters  A,B  in  the  feed- 
back path,  will  be  derived. 

1.    The  System's  Block  Diagram 

Assume  that  the  system  under  consideration  has  a  block  diagram 
as  in  figure  3-1. 
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It  turns  out  that  the  transfer  function  C/R  is: 

-Ts 

C  _  e        N(s) 

R       D(s)   +  e"Ts(As  +  B)N(s)  (IH-1) 

Therefore,    let  As  +  B  =   AN, (s)   +  BN?(s)    and   assume   that   the   coef f iciencts 
of  D(s)    are   [dklk=0    >    °f  N2(s)    are   [cklk=0    ,    and   of  N^s)    are  Cbk3k=0      » 
where   n   is    the   order   of  D(s)    and  I    is    the  order   of   N(s),    increased   by 
one,    as   N    (s)   =   s    •    N(s). 

2 .  Characteristic    Equation 

It    follows    from  the   above  discussion  that   the   system's   charact- 
eristic   equation  can  be  written   in   the    following    form: 
n  I 

k=0  k=0 

where  a,     =  Ab     +  Be  (III- 3) 

3.  Derivation  of  a  Set  of  Equations  That  Can  Provide  Constant  "x  " 
or  Constant  "g"  Curves  in  the  "A,B"  Plane.  n 

Starting  now  from  equation  (III-2)  one  can  express  the  complex 

quantity  "s"  in  terms  of  (X   and  £  and  then  attempt  to  get  two  equations 

of  the  form: 

A'B1(UJn,C)  +  B-C1(xn^)  +  D10JJn,C)  -  0  (III-4) 

A-B2(uJn,C)  +  B.C2(Xn,^)  +  D2C.Un,C)  =  0 

One  can  see  that  the  solution  follows  Siljak's  basic  steps.   In  this  case, 
then,  a  solution  for  A  and  B  is  possible,  for  a  given  (x  ,£)  pair.   Easily 
then  the  constant  "x  "  and  the  constant  "£"  curves  will  follow  by  keeping 
one  of  the  two  variables  constant,  and  varying  the  other  over  a  desired 
range  of  values. 

Shown  in  Appendix  B  is  a  way  a  deriving  equations  in  the  form 
of  the  ones  in  (III-4) .   In  fact  one  gets: 
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k=°  k=°  k=°        (Ill-Sa) 

I 

4*^  IVV^^'V  +  Vl«'V)]  + 


k=0 

I 


■[•"^Iv'J'O^k.v  +  WC'.VX 


+ 


(III-5b) 


k=0 


[ZVn(^\«)+  W^)]   =   ° 


k=0 


where  a    ,  3  ,  c,  are  constant  coefficients  of  the  characteristic  equation 
of  the  system  as  mentioned  before,  and  T,     .  ,  >l    are  given  from  the 
following  recurrence  relations: 


and  subject  to  the  following  initial  conditions 


x-i =  -  Ul<2  siiT-"Jl<2]  +z  coi^J1-^] 


x     =  cos^  /i-r2l 

o  L      riV       '  J 


xi  =    Ul<2  s±iTdjJl<2]  -  c  co{^Jl<2]) 


(ITI-6) 


'-1   '   < 


T        =      1 
o 


*i   -  -c 


By   comparing   equations    (III-4)    and    (III-5), 

i 

B.      =   eTJ^  7b.  CO  \  <£,UL    ) 

1  Zj  k   n     k  s      n 

k=0 


(III-7) 


A3 


and  the  rest  of  the  values  for  [B.,C.,D.].  ,  „  follow  easily, 

1      1      1   i=l,2 

B-C-   -   B  C     ^   0,    from  equation    (III-4)    follows: 


So   if 


A  =    (C1D2   -    C^)   v    (BlC2 


*2C1> 


(B2D1    -    D2B1)   v    (B1C2    -    B2CL) 


(III-8) 


4.    Derivation  of  an  Equation  That  Can  Provide  Constant  "Q""  Curves. 
The  real  roots  of  equation  (III-2)  may  be  found  by  substituting 


s  =  o" ,  the  required  root.   It  follows  then: 
n       I 

k    "  -TO"   //  _k 
\ 
k=0     k=0 


L\a   +Le 


where 


or 


a.  =  A*b.  +  B«c, 
k      k      k 


n       I  I 

V\   k   \"  -To  „   k   V  -TO     k 

L\°  +Le      A\°  +Le      BckCT 

k=0  k=0  k=0 


(IH-3) 


-To 
and    after   bringing  A,B,    and    e  outside   of   the    sumniation   operator,   we   get 


finally: 


-To  \",       k~l 

,e      L\°  J  + 

k=0 


Z  n 

k"1     ,   \\      k 

k=0  k=0 


(III-9) 


5.    Computer  Program 

In  Appendix  D,  a  computer  program  able  to  draw  constant  o,  £ 
and  0J   curves  or  combinations  of  them  is  listed.   Computations  take  place 
as  follows  (Notice  that  reference  is  made  to  equations  of  Appendix  B): 

a.    The  Case  of  Constant  "£"  and  "uu  "  Curves 

*        n 

(1)   B  ,  C  ,  D  ,  B  ,  C?,  and  D-  are  calculated  from  Equa- 
tions (B-15)  for  the  constant  quantity  of  the  curve  wanted.   For  each 
value  of  the  running  variable,  UU   ,  X  ,  X   ,  T  and  T    are  calculated 
by  iteration  from  the  appropriate  recursion  and  initial  value  formulas 
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as  mentioned  before  and  multiplied  by  the  proper  b  ,  c   and  d.  as  necessary. 
The  quantity   e   n^  is  also  calculated  to  form  [B.,C.]      . 

(2)   By  using  equations  (b-16)  a  point  P(A,B)  is  drawn  on 
the  A,B  plane  for  the  constant  value  curve  in  discussion,  whenever  the 
running  variable  takes  one  of  the  values  we  specify  it  to  span. 

b.  The  Case  of  Constant  "a"  Curves 

As  one  can  see,  equation  (III-9)  represents  a  straight 

line  in  the  A,B  plane  for  given  "a".   So  in  the  computer  program  of 

Appendix  D,  the  quantities, 

I  i  n 

-to    r,     k  \   k   ,  Y,     k 

e     >  l\°  >  L\°    and  L\° 

k=0    k=0       k=0 
are  calculated  for  selected  "a",  and  constant  o    lines  are  plotted  by 
using  appropriate  intercepts. 

c.  Basic  Inputs 

In  addition  to  the  types  of  curves,  the  selected  values 
needed,  the  span  of  the  running  variable  and  some  other  information  nec- 
essary to  the  drawing  unit,  one  must  input  the  order  of  the  character- 
istic equation  (which  coincides  with  the  order  of  D(s)  for  the  most 

usual  case  when  order  D(s)  >  N(s)),  the  maximum  degree  of  the  exponent 

-Ts 
in  "s"  that  the  time  delay  term,  (e    ),  appears  in  the  characteristic 

equation,  the  time  delay  constant  "T"  and  the  coefficients  of  the  system 

b.  ,  c,  ,  and  d.  . 

k'   k       k 

B.    OTHER  PRACTICAL  CONFIGURATIONS  -  GENERALIZATION 

1.    Other  Possible  Time-Delay  and  Plant  Configurations  for  Feed- 
Back  Compensation  (A  and  B  Shown  in  the  Feed-Back  Loop) 

In  order  to  get  a  generalization  of  the  parameter  plane  equa- 
tions for  a  system  with  time  delay,  one  should  devise  possible  configurations 
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of  the  plant,  the  compensation  unit  and  the  time  delay  and  after  solving 
for  the  parameter  plane  equations  for  each  case,  should  then  try  to  for- 
see  how  these  equations  would  look  for  the  general  system  with  time  delay. 

a.  In  Figure  3-2  the  first  of  a  set  of  three  configurations 
is  shown.  But  noting  that  the  loop  gain  of  this  system  is  exactly  the 
same  as  for  the  system  discussed  in  section  IIIA,  one  shouldn't  expect 

any  differences  in  the  characteristic  equations  of  these  two  systems. 

-Ts 
Therefore,  the  characteristic  equation  is  D(s)  +  N(s)(As  +  B)e    =  0. 

b.  The  next  system  to  consider  is  the  one  shown  in  Figure  3-3. 

Since  for  this  system: 

N(s) 
C  =  D(s) 

■i_  i+!8  ["•-*  +  "] 

it    follows   that   the  characteristic   equation    for   this    system  will  be: 

-Ts 
D(s)   +   B   N(s)   +   e        As    N(s)    -   0 

So  the  basic  difference  with  the  original  system  of  section  IIIA,  exists 

in  the  fact  that  the  term  multiplying  the  variable  B  is  N(s)  for  the 

-Ts 
latter  case  instead  of  N(s).e     for  the  former.   As  a  consequence  of 

that  and  going  over  the  procedure  of  section  IIIA,  one  can  see  that  no 

e   n^  term  appears  in  the  set  of  equations  (III-5)  as  far  as  C,  and  C~ 

are  concerned.   Furthermore,  only  T.  ,U.  terms  will  show  up  in  the  solution 

for  C,  and  C~„   As  far  as  B  ,  B  ,  and  D  ,  D~  are  concerned,  there  will  be 

no  change  from  the  solution  of  section  IIIA. 

c.  The  last  system  to  be  considered  for  the  feed-back  com- 
pensation case  with  A  and  B  in  the  feed-back  loop,  is  the  one  shown  in 
Figure  3-4.    In  this  case: 

C Ms) 


R  -Ts 

D(s)  +  N(s)(As  +  Be    ) 
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Therefore,  the  characteristic  equation  is: 

D(s)  +  AsN(s)  +  e~TSBN(s)  =  0 
Again  by  comparison  with  the  characteristic  equation  of  the  system  dis- 
cussed in  Section  IIIA,  one  can  see  that  the  difference  this  time  exists 

in  the  term  multiplying  the  "A"  variable  parameter.   In  fact,  instead  of 

-Ts 
A(s*N(s)*e   )  it  is  A(s'N(s)).   So  following  again  the  procedure  in 

-TJJ  C 
Section  IIIA,  one  will  get  no  e    n   term  appearing  in  the  set  of  equa- 
tions (III-5)  as  far  as  B,  and  B_  are  concerned  this  time.   Furthermore, 

only  T,  ,  U,  terms  will  show  up  in  the  solution  for  B..  and  B_.    The  solu- 
te  k  12 

tions  for  C1 ,  C  ,  and  D1 ,  D_  will  remain  the  same  as  for  the  system  of 

Section  IIIA. 

d.    Conclusion 

The  solutions  for  the  parameter  plane  equations  for  the 

systems  described  in  paragraphs  11131b,  IIIBlc,  remained  essentially  the 

same  as  for  the  system  in  section  IIIA,  with  the  only  difference  of  minor 

changes  on  the  C.  or  B.  coefficients  (i  -    1.2),  whenever  the  time  delay 

factor  doesn't  show  in  the  characteristic  equation  to  multiply  one  of 

the  two  variable  parameters  A,  B.   These  results  are  tabulated  in  Table 

3  -  1. 

2 .    Possible  Time  Delay  and  Plant  Configurations  for  the  Case  of  a 
Lead-Lag  Compensation  on  a  Fixed  Feed-Back  System. 

a.    Consider  first  the  system  as  in  Figure  3-5.   The  transfer 

function  for  this  system  is: 

-Ts 


(s  +  A)N(s) 
C  =  (s  +  B)D(s) 
R  "'  T    (s  +  A)  N(s)   -Ts 

1  +kCs  +  ^D(s)'e 


(s+A)'N(s)-e"Ts 


sD(s)  +  e"Tsk-s'N(s)  +  BD(s)  +  e"TsA(K»N(s) ) 
k         k  k 


(111-10) 
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By  examining  closer   the  characteristic  equation  of  this  system  the 
following  facts  should  be  noted: 

(1)  As  far  as  the  terms  multiplying  the  A  and  B  variable 
parameters  are  concerned  the  situation  is  similar  to  the  one  examined  in 
paragraph  IIIBlb  with  the  system  shown  in  Figure  3-3. 

(2)  For  the  first  time,  terms  independent  of  A  and  B  and 

-Ts  tt 

of  the  form  "e    ^(polynomial  in  s)  are  encountered. 

b.  In  the  case  where  the  time  delay  appears  in  the  feed-back 
path  with  all  the  other  components  of  the  system  being  as  shown  in  Figure 
3-5,  the  loop  gain,  and  so  the  characteristic  equation,  remain  unchanged. 

c.  The  characteristic  equation  of  the  system  in  Figure  3-5 

can  be  formulated  in  a  more  elegant  way  as  follows: 

n       n  n         n 

\\   k   V  -Ts,/  k   _  V   k     V  -Ts,   k   n  /TTT11N 

}AS  +lc    dks  +BAcks  +Ab    V  =  0       (III-U) 

k=0     k=0  k=0       k=0 

where  d  ,  d  ,  c  ,  and  b   are  the  constant  terms  encountered  in  the  follow- 
k   k   k       k 

ing  polynomials  in  "s",  respectively: 

sD(s),  ksN(s),  D(s),  and  kN(s). 
The  fact  that  the  summation  operators  shown  in  the  equation  (III-ll)  are 
all  operating  in  a  range  from  "0"  to  "n"  is  perfectly  legitimate  and  what 
is  meant  by  that  is  illustrated  as  follows  and  in  relation  to  the  system 
of  Figure  3-5. 

(1)  Consider:   k  =1  and  N(s)        2s  +  3 

D(s)     .3,2    ,      ,_ 

v  '    4s  +  5s  +  6s  +7 

where  order  D(s)  >  order  N(s),  as  usually  happens  with  most  plants. 

(2)  It  follows  that: 

(a)  sD(s)  =  4s  +  5s  +  6s  +  7s  +  0s 

(b)  ksN(s)=  2s2  +  3s1  +  0s° 
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3     2     10 

(c)  D(s)   =  4s  +  5s  +  6s  +  7s 

(d)  kN(s)  =  2s1  +  3s? 

where  polynomial  of  max(order)  =  4,  is  sD(s).   But  for  the  purpose  of 

neater  summation  operations  in  a  digital  computer  one  could  form  the 

characteristic  equation  as  implied  by  equation  (III-ll),  as  follows: 

n=4     n=4  n=4       n=4 

)d,  s  +  )e   d,'s   +Byc,s  +  A  )e   b,  s   =0     (III- 11a) 
Lk     Z_i     k       Z_jk       L>  k 

k=0     k=0  k=0       k=0 

where  d  ,  d'   c,  ,  b,  for  the  system  are  tabulated  in  Table  3-2,  and  n=4 
k   k   k   k 

is  implied  by  the  fact  that  the  polynomial  of  maximum  order  for  this  case  !.'. 
of  order  4. 

TABLE  3-2 

Polynomial  .      k 

4       3       2       10 

sD(s)      d.  =  4   d0  =  5   d0-6   d..  =  7   d   =0 
4       3       2       1       o 

k  sN(s)     d!  =  0  d!  =  0  d'  -  2  d'  =  3  d'  =  0 

4       3       2       1       o 

D(s)     c.  =  0   c„  =  4  c_  =  5   c,  =  6   c   =7 

4       3       2       1       o 

kN(s)     b.-0  b_  =  0  b_  =  0  bn=2   b   -3 
4       3       2       1       o 

3 .    Generalization 

a„    Having  in  mind  the  form  of  equation  (III-ll)  as  well  as 

the  results  tabulated  in  Table  3-1,  one  can  deduce  that  a  more  general 

form  can  be  given  to  the  characteristic  equations  of  all  systems  discussed 

up  to  now  in  section  III  and  this  is  as  follows: 

n       n  n  n 


L\s   +  Le  dks   +  B  b  V   +  A  Le  V   =  ° 

k=0     k=0  k=0  k=0        ClII-12) 

where  d  ,  d' ,  c,  ,  and  b   are  defined  as  in  the  previous  illustrative  ex- 

K     K.     K  K. 

ample  of  this  section  (p.   ),  and  6,  ,6  ,6,  are  variables  that  take  on 

d2    c    b 

values  of  0  or  1  depending  on  the  system  under  examination. 
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b.  Constant  "a"  Curves 

For  s  =  6,  one  gets  by  direct  substitution  from  (111-12) 
the  following  equation  representing  as  usual,  a  family  of  straight  lines: 

A  [£e-6b-T'abkok]  +  B  [2e-6c-T-acko-k]  +  [V(dk+  e"6^  'T*  °d  ■  )ak]  =  0  (111-13) 
k=0  k=0  k=0 

c.  Constant  T"  or  Constant  "jj  "  Curves 

*  n 

It  has  already  been  defined  that: 
0 


i(i=d2,c,b)    I  l 

If  now,  6'  is  defined  as: 
J 

6j(j=d2,c,b)={°r  • 

it  follows  that: 

6.  +  6!  s  1  (III- 14) 

1    J 

So  by  using  6.  and  6!  one  can  selectively  maintain  or  discard  the  differ- 
ent terms  that  must  or  must  not  show  up  in  the  parameter  plane  solution 
of  a  particular  system  if  this  system  has  a  characteristic  equation  of 
the  form  as  in  (111-12). 

(1)  Example.  By  examining  the  solutions  for  B.  and  B~ 
in  table  3-1,  it  follows  that  for  the  systems  discussed  up  to  now,  B. 
takes  one  of  the  forms:       . 

(a)  Bi  - 1\\\ 

or    (b)      BL   =   e      vS  Zjbkxn  xR 

k=0 

This  statement  can  be  reduced  into  one  single  equation  by  using  (111-14), 

and  the  fact  that  B.,  is  related  to  b,  's  and  therefore  to  the  equation 

1  k 

5,-f6*  =1    (observe  6,  to  be  related  with  "b,"  coefficients  in  (111-12). 
d   d  b  k 
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(c)      B.    =  e"6b,T,^*C.6    .  yb  M  \     +&  .yb  (JJ  kT      (in-15a) 
1  bZ_iknkb/_knk 


k=0  k=0 


(d)  Similarly: 

k  =  0 
n  K   U  (III-15b) 

k=0 
Solutions  for  C,  ,  C_  and  D,  ,  D9  follow  similar  patterns. 

(e)  It  is  of  interest  to  note  that  if  6.  =  1  (6!=0), 
the  presence  of  X  's  and  X,  , ' s  is  implied  in  the  relevant  expression 
for  B  ,  B  ,  C  ,  C  ,  or  D  ,  D  .   If  in  the  contrary,  6 !  =  1  (6.=0),  the 
presence  of T,  's  and  T,   's  is  implied. 

(2)   Derivation  of  Equations 

For  one  to  actually  solve  for  the  parameter  plane 
equations,  one  must  make  use  of  the  fact  (as  in  section  IIIA),  that: 

s  =  'JJ  (cos0  +  jsin9)  =  JJ  eJ 
n  n 

and  apply  it  in  equation  (111-12)  to  get: 
n  n 

\j   k/    ,  a    •   •   ,an   Vji/  -OJo'T*  (cos9  +  jsin9)a:ri)  k   jk9 
)d,U)   (cos  k9  +  i  sin  k9 )  +  )d,'(e   z  J      /JJ   •eJ   + 

Li  k  n  J  L>  ky  n 

k=0  k=0 

n 


'/_.   -6cT(cos9    +  jsin9)iJnN>  k   jk9 


k=0 
n 

^Ae-6bT(cos9  +  jsin9)V^keJk9  =  0 

k=0 
which,  with  appropriate  manipulation  inside  the  summation  operators  in- 
volving d ' ,  c   and  b  ,  yields  the  result  shown  in  equation  (III- 16). 

K.     iC         K. 

(See  next  page.) 
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)d,UO        (cos  k9   +  j    sin  k9  )   + 

Lk  n 


k=0 


6d2,T,(VC  Yd^JUnk  {[V    •    cos(k9    +  uun-T-C)   +  6^    'cosCke)]   + 
k=0 

j  |_6      •    sin(k9    +  LUn«T'C)   +  6^2-sin(k9)]j      + 


B^c'T'i'n'C    •      w  k  |[  6      .    cos(k9    +  w    .T.^)   +51    .Cos(k9)]   + 
L  k   ii      iL    c  n        ^  c  J 

j  I  6      •    sin(k9    +  uu    •*!>£)   +  61    •sinfcS)]!      + 


k=0 


n 


Ae6b'T-Xn'C  jbkuunk  {{Y     •    cos(k9   +  m   -T'O   +  6b    .cos(k9)'j  + 
k=0 

j  I  6      •    sin(k9   +  'X'    »T»C)   +  &'    •sin(k9)Jj 


(111-16) 


From  this  stage  by  using  appropriate  substitutions  for  cosines  and  sines 
with  arguments  (k9 )  or  (kO  +  UU  •T*£),  in  terms  ofX  's,Y  's,  T  's,  and 
U,  '  s  and  their  recurrence  relations,  exactly  as  defined   derived  and 
used  in  section  IIIA,  one  can  form  the  following  equations  after  separ- 
ating reals  from  imaginaries: 

n  n 

ydw  kT  +  e6d2'T,'!Jn,C  d^   k(6   .x   +  6  '  T  )   + 
L>  k  n  k  Li  k  n  v  d2  k    d2  k 

k=0  k=0 

n 

Beuc  i  an  C  )c  w   K(6  y  +  6»T,  )  + 
Li  k  n   c  k    ck 

k=0 

n 

k=0 


6b.T-JJn-C  yb  a,  k(5X   +  6'T  )  =  0  (III-17a) 

L,  k  n   b  K    b  k 


(see  next  page) 
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and: 


Z_iknk  Z_ikn        cl2   k  d2   k 


Be 


6r»T*uun«£  k,c-  c,„   , 

n  *»   /  c,  tu        0   V,    +0  U.)   + 
Z-j  k   n        c   k  c  k 


Ae6b,T*,JVC  Vb  u,  k(6   v     +  6'U,  )    =   0 

L  k  n       b  k  b  k 


(III-17b) 


Equation    (III-17b)    can  be   expressed    in   terms    of         's    and    T    '  s    instead    of 

K.  K. 

U,  's  and   ,  's,  in  exactly  the  same  manner  as  in  section  IIIA  by  using 

K  K. 

equations  (B-9c)  and  (B-12c)  of  Appendix  B.   So  the  second  equation  of 

the  parameter  plane  system  becomes: 
n  n 


L 
k=0 


A'>Tk  +  Tk-1>    +  c6d2'T"'Xn'CIdk'nk[6d,«\  +   Vl»  V*  +  Tk-1>]   + 


k=0 


Be 


6c-T^n.CVr        1 


k=0 


Ae 


6b'T'ln'c  LWl^h  +  W  +  ab«Tk  +  Tk-!>]  ■  ° 


(III-17c) 


k=0 


For  easier  reference,  recurrence  equation  as  well  as  special  symbol  re- 
lations are  listed  below: 

.   (a)   Tk+1  -  -2{Tk  -  tl 

(b)  Jk¥l   =  -2ZXk   -    Xk_x 

(c)  Initial  condition  relations  usually  necessary 
in  digital  computer  calculations: 

(i)  T_x   =  -C,  TQ  -  1,  TL  =  -C 


(ii)  x_1  =- 


^/l-C^inCTx^l-C^+CcosCTJjJl-C2)] 


x     =  cos(ti  Vi-r2) 

o  nv   ^ 
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Xl   =Jl<2sin(TJJJl<2>    -  CcosCTJ^l-C2) 


(iii)     6!  +6.   =  1 
J       1 

i    =    J    =    (d2,c,b) 

(3)  Parameter  Plane  Curves.   Constant  "JJ  "  and/or  con- 
n 

stant  "£"  curves  can  be  drawn  in  the  parameter  plane  in  exactly  the  same 
manner  as  discussed  in  section  I-B  (Siljak's  approach). 

(4)  Conclusions . 

(a)  By  comparison  of  equations   (III- 17)  and  the  sol- 
utions for  B. ,C1 ,D. ,B_ ,C_ ,D_  tabulated  in  Table  3-1,  one  can  observe  that 
the  derived  equations  (111-17)  completely  cover  each  case  and  furthermore 
they  cover  the  case  of  the  system  shown  in  Figure  3-5. 

(b)  The  equations  (111-17),  are  valid  for  the  case 
that  the  terms  encountered  in  the  characteristic  equations  examined,  and 
involving  the  variable  parameters  A  and  B  are  either  of  the  form: 


-Ts 
(Ae    )• (Polynomial. (s))+  B* (Polynomials (s)) 


or  of  the  form: 

-Ts 
A* (Polynomial- (s) )  +  (Be    )• (Polynomial, (s) ) 

This  was  because  of  the  specific  features  of  the  systems  examined.   But 
since  the  terms  not  involving  A  and  B  in  equation  (III-ll)  are  of  the 

form: 

-Ts 
e    • (Polynomial- (s) )  +  Polynomial, (s)  , 

the  possibility  exists  for  a  more  general  case  with  A,B  as  well  as  AB 
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appearing  in  the  characteristic  equation: 

F(s)  +  F'(s)e"Ts  +  J"b(s)  +  B'(s)e"Ts]+  b[c(s)  +  C  (s)e"Ts]+ 

AbTd(s)  +  D,(s)e"Tsl  =  0  (111-18) 

The  parameter  plane  analysis  of  such  a  general  system  follows  next. 

C.    NON- LINEAR  VARIABLE  PARAMETER  PRODUCTS  INCLUDED  IN  THE  CHARACTER- 
ISTIC EQUATION 

1.  In  figure  3-6  a  two  variable  parameter  system  is  shown,  which 
will  result  in  a  characteristic  equation,  containing  the  AB  product  of 
these  two  parameters,  as  follows: 

F(s)  =  [d(s)  +  sN(s)e"Ts]  +  AjN(s)e"Ts]  +  b[  s2N(s)e"Ts]  +  AB[sN(s)e"Ts] 

=  0 

or, 

n+1  n  n+2  n+1 

F(s)  =^(fksk  +  f^ske"Ts)  +  A  J^sV  S)  +  B^(c^ske"TS)  +  AB  ^(dksVTs) 

k=0  k=0  k-0  k=0 

=  0  (111-19) 

2.  However,  in  order  to  gain  more  generality  so  as  to  apply  the 

results  into  a  broader  system  classification,  it  is  wise  to  solve  for  the 

parameter  plane  equations  of  the  following  characteristic  equation: 
n  I 

k=0  k=0 

m  p 

^B(ck  +  c^e"Ts)sk  +^AB(dk  +  dke~Ts)Sk  =  0         (111-20) 
k=0  k=0 

and  by  arbitrarily  choosing  n  =  max(n,  i ,  m,  p)  as  well  as  by  filling 

non-existent  coefficients  of  lower  degree  polynomials  by  zeros: 
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n  n 


k=0  .       k=0 

n  n 

bY(c.    +  c/e"Ts)sk  +  AB  }(d,    +  d/e"Ts)sk     =     0 
Li     k  k  L,     k  k 

k=0  k=0 

3.         The   solution    for   constant  UU      and   t   curves    is    again  based    on   the 

n 


equations  : 

s   =  UU    (cos  9   +  i    sin  9)   =  au  e-1 

n  J  n 

and  x=   cos(k9    +  UU  T£),      y   =    sin(k9    +  UU   T£) 

T  =   cos   k9  ,      U  =    sin  k9 

K.  K. 

(escorted   by   their   recursion   relations),    being      substituted    in  equation 

(111-20),    and    thereafter  on   a    separation   of   reals    and    imaginaries. 


4.         Actually,    the    set   of    the    two   equations    is: 

(1-11) 


AB     +  BC      +  ABD      =    -F 


AB„   +  BC2   +  ABD      -    -F 


where: 

n 


B     =   >b,u:     T,    +  e     nS>   /b'uu      X 
1       L,  k   a     k  Li  k   n     k 


k=0  k=0 

n  n 

C.  =  )c   uu  T     +   eTln£   yc/uu  kX 
1       /_.  k   n  k                      Z_iknk 

k=0  k=0 

n  n 

D.  =  )d,UJ  kT,     +  eTJJnC  ydiu,  kx 
1       Li  k   n  k                     ^_.knk 

k=0  k=0 

n  n 

Vr      k™  Tx-nr   ~       k 

Fl  =  LVn  \  +  e           L"k\   \ 

k=0  k=0 
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n  n 

B,   =  Yb.au  k(CT.  +  T.     .)   +  e^   W'ju  k(^1    +  X      . ) 

2       Z_i  k   n     s    k  k-1                        t-,  k   n     s  k          k-1/ 

k=0  k=0 

n  n 


k=0 
n 


k=0 
n 


D2   "  ZVn^k  +  Tk-1>  +  eTV  Adk^nk«XR  +Xk-1> 
k=0  k=0 

n  n 

F2  "  ZA'U>Tk  +  Tk-1>  +  ^  Z£klnk("k  +  Vl> 
k=0  k=0 

5.    According  to  Hollister,  [2],  then: 


B  = 


(A   +  A  )  ±  /(A   +  A  )   -  4A  A 
^HFD    CB'   V^TD    CB'      CD  FB 


2A, 


CD 


and  after  B  is  defined 


A  = 


F1  +  C  B   F2  +  C2B 

B   +  D  B  "  B„  +  D  B 


where : 


A 


CD 


Cl    C2 

Dl    D2 


**I 


Dl     D2 


,  etc 
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IV.   STABILITY  ANALYSIS  IN  THE  PARAMETER  PLANE 

In  this  section  the  stability  concepts  in  the  parameter  plane  are 
built  very  qualitatively,  at  first  for  a  delay-free  system.   These  con- 
cepts are  based  on  ideas  similar  to  the  ones  used  in  conformal  mapping. 

For  the  case  of  systems  with  time  delay,  Karmarkar's  method,  based 
on  Cauchy's  principle  of  Argument  or  Michailov's  test  and  its  extension 
by  Pontryagin,  is  followed  and  used  to  examine  the  stability  region  of 
the  system  shown  in  figure  6-1. 

A.    SYSTEMS  WITH  NO  TRANSPORT  LAG 

1.  Given  a  certain  polynomial  CE(s)  in  "s",  and  of  nth  degree.   In 

general  out  of  the  n  roots  of  this  polynomial , r  will  lie  in  the  shaded 

— — —  J  /     s 

region  of  figure  4-la  and  r   in  the  unshaded  one,  but 

n  =  r  +  r  (IV-1) 

u    s 

So  the  entire  "s"  space  is  divided  into  two  sub-spaces,  by  the  s  =  ji1 
line  and  the  o  =  0  point  which  in  a  sense  are  their  complex  and  real  root 
boundaries . 

2.  Now  considering  the  closed  contour  c,  /consisting  of  c.  and 

c_  (-  the  j(JU  axis)  as  shown  in  figure  4-la},  we  have   to  realize  that  it 

maps  in  the  AB  space  by  appropriate  mapping  of  its  boundary.   But  since 

the  CE(s)  polynomial  serves  as  the  mapping  function,  it  follows  that  for 

different  characteristic  equations  the  O   =  0   and  s  =  jx  boundaries  will 

map  differently.   It  should  also  be  noted  that  these  boundaries  will 

generally  decompose  the  whole  AB  space  into  a  finite  number  6  of  regions 

R.  ,  each  of  which  will  contain  i  roots  in  it.   So: 
6 

) i  =  n  •       (IV-2) 

1 
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3.  The  whole  purpose  of  the  stability  analysis  is  to  define  the 
number  of  roots  existing  in  each  of  the  regions  R. .  (In  other  words,  "to 
determine  the  region  index'.') 

4 .  Contour  Shading  (The  Shading  Rule) . 

Consider  the  mapping  process  of  figure  4-lc,  where  the  complex 
root  boundary  is  mapped  in  the  AB  plane.   So  if  the  mapping  device  lies 
at  point  M,  it  follows  that  CE(s)  has  a  root  s  =  jut)  .   In  other  words, 
"s"  maps  into  "M"  via  the  mapping  function  CE(s)  =  0. 

Considering  now  two  points  in  the  M  neighborhood,  which  lie  in 
either  side  of  the  £  =  0  boundary,  one  realizes  that  while  moving  from 
point  M,  to  M„ ,  at  least  one  complex  root   in  the  s  plane,  moves  from 
s.  to  s„  or  vice  versa.   It  is  therefore  necessary  to  orient  the  boundary 
regions,  according  to  the  following  rules: 

a.  In  the  s  plane,  a  right  handed  rectangular  system  (t,n) 
is  established,  where  t  is  the  positive  tangent  direction  and  n  is  the 
normal  to  it  in  the  appropriate  direction.   The  boundary,  then,  is  shaded 
always  in  the  direction  of  positive  n.   (See  figure  4-lc,  upper). 

Since  the  complex  root  boundary  in  the  s-plane  is  s  =  j JJ , 
positive  t  is  determined  by  direction  of  increase  of  the  variable  working 
along  it. 

b.  In  the  AB  plane,  although  the  (t,n)  system  is  always  right- 
handed  in  the  s-plane,  this  is  not  the  case  for  the  AB  plane.   So  depend- 
ing on  the  particular  system  under  examination  the  positive  n  direction 
can  be  either  side  of  t.   The  positive  t  direction  is  determined  as  before, 
So  for: 

(1)   (t,n)  Right-Handed;   The  positive  "n"  direction  is  to 
the  left  of  positive  "t",  and  the  boundary  is  shaded  to  the  left  of  pos- 
itive "t". 
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(2)   (t,n)  Lef t-Handed :   Since  positive  "n"  will  be  direct- 
ed to  the  right  of  positive  "t",  the  boundary  is  shaded  to  the  right  of 
positive  "t". 

c.  The  Significance  of  Shading. 

The  shaded  portions  of  the  regions  in  both  the  s  and  AB 
domains  correspond,  in  the  specific  case  examined,  to  the  side  of  the 
domain  where  all  the  complex  roots  exist.   That  is  to  the  stable  region 
since  the  boundary  in  "s"  was  the  absolute  stability  boundary. 

d.  Rule  to  Define  the  (t,n)  Co-ordinate  System's  Orientation 

(1)  Consider  the  polynomial 

F(s)  =  CE(s)  (IV-3a) 

(2)  Substitute  s  =  O   +   j'JU,  to  get 

F(a,uu)  =  R(a,yj)  +jl(a,uu)  (IV-3b) 

(3)  After    forming    the   Jacobian 

J  =  J  [(R,I)/<A.B)]  =  (||)(|i)  -  (|fX|i)  (IV-3C) 

decide  as  follows: 

(a)  for  J  >  0,  (t,n)  is  right-handed,  therefore 
shade  to  the  left  of  positive  "t". 

(b)  for  J  <  0,  (t,n)  is  left-handed,  therefore 
shade  to  the  right  of  positive  "t". 

e.  Similar  rules  apply  for  the  real  root  boundary,  (the 

0"  =  0  curve  for  this  case).   However,  the  following  rule  applies  success- 
fully: 

(1)   Since  at  the  point  of  intersection  of  s  =  0"  with  the 
s  =  juo  curve  uu  changes  sign,  so  does  the  Jacobian  J.   Therefore,  in  the 
AB  plane  the  0"  =  0  curve  must  change  direction  of  shading  at  the  point  of 
intersection  with  the  £  =  0  curve. 
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(2)   In  the  stability  region  neighborhood  of  the  image  0 
in  the  s-domain,  the  shading  of  the  o   =  0   and  £  =  0  curves  must  coincide, 
since  this  is  the  case  for  the  s-domain  shading  also. 

5.    The  stability  region  in  the  AB  plane,  for  a  delay- free  system, 
has  now  been  established.   For  more  details  one  could  refer  to  reference 

Ei]. 


B.    SYSTEMS  WITH  TIME  DELAY 

1.  In  this  subsection,  Karmarkar's  approach,  L6],  to  define  if  the 
domain  of  maximum  number  of  roots  is  also  the  stability  domain,  will  be 
followed  since  it  provides  in  a  single  plot  the  stability  domain  as  well 
as  the  test. 

2.  Karmarkar's  stability  theorem  states  the  following: 
"Consider  CE(s)  to  be  a  transcendental  polynomial: 


m,n 


CE(s)  =  )  a   sPerS 
Li     pr 


P,r 


where: 


a   are  real  coefficients 
pr 

a   4-   0  (the  principal  term) 
mn  r      r 

p,n  are  intergers,  and 


I 


p  rs     0   1 
a   s  e   =   s   s 
pr        L 


...s«] 


(IV-4) 


aoo    aoi  •"  a 
aio    an  *••  a 


a  /-i  a  0  ...  a 

mO   ml 


m 

r  osi 

e 

On 

Is 
e 

in 

mn_ 

'ns 
e 

(IV-5) 


For  CE(s)  to  be  stable  it  is  necessary  and  sufficient  to  find  an  interger 

K  so  that,  starting  with  the  line  B  =  b  ,  the  total  number  "X"  of  alter- 

o 

nating  intersection  of  the  T  =  0   curve  with  the  A  =  a   and  B  =  b   lines 

v  o  o 
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is: 

X   =   4kn  +  m  (IV-6a) 

when  uu  ranges  over  the  interval 

0  £  uu  £  2nk  +  e/2,  for  0  =£  e  £  n/2, 
or  for  s=tt6,   0  £  uu  ^  2iT(k  +  6/4)  for  0  £  6  £1.  (IV-6b) 

In  figure  (4-2),  the  results  of  this  analysis  for  the  system 

CE(s)  = (s  +  2s   +  s)eS  +  as  +  b  =  0 
is  shown.   For  this  particular  system  one  has  m  =  3  and  n  =  1.   So  since 
"uo"  was  varied  in  the  range  uu  =  (small  v  14)  and  therefore  : 

uu  =  14  <  2n(k  +  1/4),   =  14.15=  uu    (where  uu  in  rad/sec) 

k=2  max 

for  k=2,  the  intersection  predicted  by  the  theorem  should  be 

\   =   [(4kn  +  m)]k=2  =  11, 

n=l 

m=3 

and  this  is  what  is  measured,  so  the  system  is  stable  in  the  domain  s 
bounded  by  the  o~  =  0  line,  which  is  the  B=0  line  as  easily  one  sees  from 
CE(s),  and  the  first  encirclement  of  the  £=0   curve  since  the  Jacobian 
(IV-3),  is  less  than  zero  for  this  case. 

3.    The  actual  derivation  of  Karmarkar's  theorem  is  based  as  already 
mentioned  in  the  following  theorems: 

a.  Mikhailov's  Criterion,  as  a  direct  application  of  the 
Cauchy  argument  principle^  7]. 

b.  Pontryagin  Theorem  on  the  zeros  of  elementary  transcenden- 
tal functions^  8] f    and  the  conditions  imposed  on  A  and  B  when  the  mapping 

* — 

operating  point  scribes  the  y  =  0  curve.  (That  is  equivalent  to  having 
roots  of  the  form  s  =  juu  in  the  s-plane).   For  details  on  his  derivation 
see  reference  [6]. 
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CONSTANT  ?  =  0  AND  <Z=0  CURVES 
FOR  A  SYSTEM  WJTH  TPELAY 


FIGURE   4-2 
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C.    APPLICATION  ON  AN  ACTUAL  SYSTEM 

1.  Consider  the  stability  of  the  system  shown  in  figure  6-1  ,  for 
different  values  of  the  time  delay  constant  T  =  0.0,  0.1,  0.2,  0.3. 

It  is  obvious  that  for  T  =  0.0  one  can  apply  methods  discussed 
in  section  IV  A,  and  for  the  rest  of  the  time  delay  constants  the  methods 
of  section  IV  B  are  necessary. 

However,  a  far  more  practical  method  of  determining  the  stability 

of  the  system  is  to  arbitrarily  pick  two  points  in  either  of  the  smallest 

*1 

two  areas  defined  by  the  closed  constant  £  =  0    curve  and  the  o   =  0 

straight  line,  and  simulate  the  system  in  "t"  domain.   If  for  one  point 
the  system  is    stable,  the  domain  it  represents  is  the  stable  region. 

2.  The  characteristic  equation  of  the  system  is: 


CE(s) 


s6  +  98.4696s5  +  4941.1680s4  +  8978.082s3  +  5086.9961s2   + 
3162.6011s    -    147.59    I   +  e"Ts  j~A(s3  +   .444s2   +   .488s)    + 
B(s2  +    .444s   +    .488)1  =    f(s)   +  e"Ts|A  g(s)   +  B  h(s)1   =   0      (IV-7) 
a.         In   order    to   obtain   constant   "a"   curves,    let    s    =  0" . 


CE(a)   =    f  (a)   +  e_TJ  [ A  g(a)   +  B  h(a)  |   =  0 

and    for    the  0=0    curve   one   gets: 

A  g(0)    +  B   h(0)    =    -f(0)    , 

or        A  x  0  +  B   x    (.488)    =    147.59 

for  finite  time  delay,  the  constant  "o  =  0'  curve  is: 

B  =  302.87  (IV-8) 


*1 

Notice  the  equivalence  of  the  £=0  and  E  =  0  curve.   In  general 

a  "E  =  a"  curve  is  the  loci  for  roots  "s  =  or  +  jlD." 
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Ts 
b.    After  multiplying  equation  (IV-7),  by  e  "  one  gets: 

Ts 


f(s)  +  e°S  I  A  g(s)  +  B  h(s)l  =  0 


(IV-7a) 


Now  since  T  is  not  an  integer,  call  T  =  tn/ tj ,  where  "tn"  and  "t<j"  are 
integers,,  and  substitute  into  (IV-7a): 

etnS,f(s')  +  e°S'  [a  g(s')  +  B  h(s')]  =   0  (IV-9a) 

(IV-9b) 


where 


s'  =   S 


td 


As    a   consequence   of   equation    (IV-5) ,    equation    (IV-9a)    can 
also   be  written   as: 

(s')O{e°S^0   x  A  +    .488b]   +  etnS  '[-147 .  59]}   + 

INif\OsT 


(s*)  |e       [ 


t    m2I  0s 
(s   )   Se 


.488A  +   .444b]   +  etflS  i  3162.601lj|   + 


,444A  +  B 


'] 


t  sT  T1 

+  e   »     I  5086.9961  If   + 
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(s')3{e°S'[A  +  0   x  Bj  +  etnS'|  8978.082 


,>4J»= 


(s')5{e°S'[o] 
(s')6{e0S'[°] 


+  etnS  |  4941.1680  if 


»1 


+  &tnS'\ 98. 4696]}   + 


+  e 


(IV-9c) 


According  to  the  definition  of  section  IV  B,  the  principal 


.      ,  m    ns. 
term,  (s  a   e   ) ,  is 
mn     ' 


6      t  s  ' 
(s1)  (1.0)e  nS 


where  m  =  6,  n  =  t   and  a    =1.0, 

n      mn 


(IV-10) 


So  for  CE(s)  to  be  stable  it  is  needed  to  find  an  integer 


k  such  that ,  for  x  *  =  (."/td)   in  the  range  0^0)'^  2nk+  -  =  tt  (2k+  t)  , 
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where: 

0  -<  e  £  ^  ,  and  0  £  6  £  l  , 

the  total  number  of  intersections,  \,  of  the  £  =  0  curve  with  the  A,B 

axes  moved  to  a  point  P(A=a  ,  B  =  b  )  is: 

o       o 

X   m   4kn  +  m  s  (4t  )k  +  6 
n 


In  Table  4-1,  the  columns  under  the  heading  "calculation", 
are  "filled  in"  according  to  these  two  formulas  giving  "X  and  ug"  for 
t  =  (1,2,3),  which  correspond  to  time  delay  constant  values  of  T  =  (.1, 
.2,  .3).   The  rest  of  the  columns  are  filled  in  relation  to  a  computer 
solution  of  the  "£=0"  curve  vSi  frequency  and  this  will  be  discussed  later 
in  this  section. 

c.    Another  fact  of  practical  importance  to  note  is  that 
karmarkar's  theorem  starts  counting  crossections  of  the  £=0  curve  with 
the  axes  through  a  testing  point  in  the  A,B  space,  from  a-0 .   But  any 
digital   machine  must  have  an  initial  value  of  uu^O . 

So  it  is  necessary  to  calculate  the  origin  of  the  £=0 
curve  for  UU  tending  to  zero.   This  limiting  point  can  be  calculated 
rather  easily  since  for  the  case  of  T=0  the  time  delay  effect  is  not 
present.   But  also  for  ID  -»  0  on  the  £=0  curve,  the  time  delay  effect  is 
also  unimportant  since: 

lim  e"Ts  =  lim  e"j^T  (IV-11) 

T-O         ttr*0 

So  the  origin  of  the  constant  £=0  curves  is  the  same  for  a  system  with  a 

variable  time-delay  constant,  and  independent  of  T. 

e.    For  the  particular  system  of  equation  (IV-7),  one  takes 

the  limit  by  operating  on  CE(s)  as  follows: 

limit  of  (£=0,  curve)    =  lim  <"E(jt) 

ul>-0    ur*0 
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But  this  operation  will  strike  out  all  the  "s"  dependent  terms  of  CE(s), 

and  therefore 

lira  CE(j'X')  =  -147.59  +  .488B 
uu-O 

and  by  equating  that  to  zero  the  solution  is  again 

B  =  302.87  (IV-12) 

One  should  anticipate  this  continuity  of  real  and  complex  root  boundary 
since  for  'JD~*0  the  time  delay  effect  is  not  present  and  therefore  the  dis- 
cussion of  section  IV  A,  is  valid. 

f.  It  remains  now  to  determine  the  value  of  the  Jacobian  in 
equation  (IV-3c).  This  Jacobian  is  in  general  a  function  of  O,  UU,  A,  B 
so  for  a  particular  o 

lira  J(ct,^,A,B)  -  lira  [  —  ~  -  —  rrj  - 

lim  — -  •  lim  — -  —  lim  — -  •  lira  — - - 
a_a  3  A   G.,a  6B   Q_a  SB   a^a  6 A 

and  because  of  the  nature  of  the  derivative,  operation,  and  the  well  be- 
haved characteristic  equation,  the  order  of  the  limit  and  derivative 
operators  can  be  interchanged.   So  for  0"-»a  =  0: 

Trn  u,  A  rn  1  5R(0,'x)   ^I(0,x)  _  aR(OA')  .  BI(0,m) 
j^u,u),a,b;  -   8A       dB       dB       sA 

Now  referring  to  equation  (IV- 3b),  R  and  I  will  be  computed 
for  T=0,  and  the  same  orientation  of  the  (t,n)  coordinate  system  will  be 
assumed  for  T  ^  0  since  for  uu  =  uu  =  0  and  a  differential  change  in  uu, 
(to  cause  uu  =  uu  +   duu  =  0  +  d'JL)) ,  the  (t,n)  orientation  shouldn't  change. 
But  once  established  the  orientation  of  the  £=0  curve  cannot  change.   So 
the  orientation  for  a  system  with  variable  time  delay  is  independent  of 
T. 

Now  taking  into  account  that  some  terms  of  R  and  I  will 
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disappear  in  the  derivative  operations,  define: 

2     2 
R   ,  (O.iD)  =  -.444x  A  -  'X  B  +  .488B 
mod 

3 

I   i(0,JJ)  =  -x  A  +  .4881A  +  ,444u)B 
mod 

to  be  regular  R,I  polynomials,  with  the  A,B  independent  terms  deleted, 

and  the  expression  for  J  is: 

J(0,x)  =  -(x)5  +  .291x3  +  .488x2  -  .238 
Considering  cases  with  x  >  0: 

lim  J(0,x)  =  -(x)5  <  0 

and  lim  J(0,x)  =  -.238  <  0 

x-0 

So  at  the  "x"  extremes  J  <  0  and  due  to  the  nature  of  the 

problem  the  same  must  occur  for  intermediate  values.   Therefore,  (t,n), 

is  Left-Handed  and  for  X  increasing  the  shading  goes  to  the  right  of  the 

£=0  curve. 

3 .    Digital  Computer  Stability  Analysis 

a.  With  all  this  information  available  the  "a  =  0"  and  "£  =  0" 
curves  are  solved  with  the  aid  of  a  digital  computer  and  the  program  of 
Appendix  D. 

b.  This  was  done  for  T  =  (0.0,  0.1,  0.2,  0.3),  and  the  re- 
suits  are  displayed  as  follows: 

(1)  Maximum  root  regions:   In  figures  from  4-3  to  4-6. 

(2)  Stability  analysis:   For  a  point  P  in  the  A,B  plane 
which  lies  in  the  general  area  of  the  system  of  section  VI  is  going  to 
to  operate, (that  is,  P  =  (A  =  5722,  B  =  10568)),  the  g  =  0  curve  was 
drawn  resulting  to  crossections  with  the  An  =  5722  and  B„  =  10568  lines 
as  shown  in  figures  4-7,4-8,4-9. 
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c.    By  using  this  information  the  rest  of  the  table  4-1  was 
filled  in  and  the  partial  conclusion  was: 

(1)  For  T  =  0.1  and  0.2,  the  system  was  found  to  be  stable, 

(2)  For  T  =  0.3,  the  "JU  range"  of  the  simulation  obtained 
showed  that  for  k=0,  k=l  the  requirement  for  stability  was  not  obtained. 

4.  General  Conclusions 

The  first  thing  to  notice  is  that  for  T  =  0.0  one  gets  the  big- 
gest maximum  root  region.   This  AB  region  is  with  no  other  reservation, 
the  stable  region  in  the  AB  plane. 

As  the  time  delay  constant  is  increased  the  maximum  root  region 
is  decreased.   But  still,  this  maximum  root  region  is  not  the  stable 
region.   According  to  Karmarkar's  criterion  the  system  with  T  =  0.3  is 
not  defined  stable  for  k=l  and  so  one  should  try  k  >  1.   These  results 
will  be  checked  out  and  evaluated  in  section  VI,  where  this  practical  de- 
sign problem  is  encountered  in  more  detail. 

5.  Notice  also  that  for  the  one  variable  parameter,  "k" ,  system  of 
figure  2-2  the  conclusion  was  that  for  increasing  time  delay  constant  the 
maximum  k  for  stability  was  decreasing.   So  in  a  sense  the  stable  "para- 
meter line"  was  decreasing. 

In  the  two  variable  parameter  system  one  has  also  a  decreasing 
area  of  stability  in  the  parameter  plane  with  increasing  T.  This  should 
be  expected  since  the  physical  situation  remains  the  same. 

As  a  result  of  that  we  should  expect  the  volume  or  the  hyper- 
volume  of  stability,  if  the  system  can  be  stable  at  all  for  a  given  time 
delay  constant   for  a  third  or  higher  order  variable  parameter  system, 
to  decrease  with  "T"  increasing. 
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TABLE  4-1 


Delay 
Constants 


n     d 


selected 


Calculation   Simulation 

A.   uu        X  uu 


max 


Remarks 


0.1 


0.2 


0.3 


1      10 


2     10 


3     10 


18   62.8(3-^) 

204.2 
rad/sec 


14   62.8(1+^) 

78.52 
rad/sec 


18   203.5 
rad/sec 


14   77.2 
rad/sec 


Stable 


Stable 


r    I 

18   62.8(l4~)!l9   77.2 
=         rad/sec 
78.52 
rad/sec 


Not  defined 
Should 
check  for 
k  >  1 
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V.   SYSTEMS  WITH  THREE  VARIABLE  PARAMETERS 

The  purpose  of  this  section  is  to  provide  a  general  way  of  analysis 
of  systems  with  three  variable  parameters  to  meet  specified  conditions. 

In  the  analysis  to  follow,  a  certain  knowledge  of  the  basis  of  de- 
scriptive geometry  is  needed,  so  for  the  reader  who  is  not  familiar  with 
the  subject,  reference  [9]  is  recommended.   However,  the  elementary  con- 
cepts needed  later  are  presented  as  well. 

The  attack  to  the  solution  of  the  most  general  case  of  a  system  with 
time  delay  is  done  indirectly.   In  fact,  a  system  with  characteristic 
equation  linear  in  A,B,r  is  first  analysed  so  that  the  concepts  of  the 
solution  are  being  built  first  for  this  simpler  case.   The  transition 
then  to  the  basic  objective  is  not  difficult. 

A.    ELEMENTS  OF  DESCRIPTIVE  GEOMETRY 

1.  Consider  a  three  dimensional  orthogonal  coordinate  system,  with 
axes  A,  B,  i  .   Shown  in  figure  5-1  is  the  set  of  points  P..  ,  P.-,,  P,-.,  P.. 

2.  The  purpose  of  the  descriptive  geometry  is  to  provide  means  of 
accurate  representation  and  measurements  of  a  three  dimensional  space  in 
a  two  dimensional  one.   So  if  for  this  purpose  the  plane  OBr  is  selected 
as  the  reference  where  the  measurements  will  take  place,  the  projections 
of  any  point  of  the  actual  system  on  the  actual  OAT  plane,  can  be  shown 
on  the  reference  without  distortion  by  just  a  90   clockwise  rotation 
around  the  Or  axis.   Figure  5-2  then  will  result  in  which  B  and  -A  as 
well  as  A  and  -B  axes  coincide.   The  projections  of  any  point  on  the  OBT 
plane  are  shown  by  "+"  and  on  the  OAT  plane  are  shown  by  "fl". 
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3 .    Points  on  the  Reference  Plane 

By  inspection  of  figures  5-1  and  5-2  one  can  note  the  following: 

a.  A  point  P  in  the  actual  coordinate  system,  can  be  repre- 
sented by  two  points  P'  and  P"  on  the  reference  plane.   So  we  have  a  one 
into  two  forward  mapping  (  and  a  two  to  one  reverse) . 

b.  These  two  points  (P',P")  that  form  the  image  of  P,  lie  on 
a  line  normal  to  the  Or  axis,  due  to  the  nature  of  the  image  generation 
process . 

c.  If  a  three  dimensional  point  lies: 
(1)  on  the  PET  plane,  its  P"  must  lie  on  the  or  axis  in  the 

descriptive  geometry  representation  (e.g.  P..  ) . 


(e.g.  P2) 


(2)  on  the  OAT  plane,  its  P1  must  lie  on  the  Or  axis  also 


(3)  on  the  A~B  or  A=-B  plane,  its  P1  and  P"  must  have  the 


same  distance  from  the  Of  axis  (e.g.  P  ) . 

(4)  on  the  Or  axis,  its  P'  and  P"  coincide  (e.g.  P/) • 

4 .    Straight  Lines  on  the  Reference  Plane 

If  now  points  P..  P  are  joined  together  to  form  the  straight 

line  they  define  in  the  actual  space,  it  is  reasonable  to  anticipate  that 

the  images  of  this  line  in  the  reference  plane  must  be  the  two  straight 

lines  that  are  defined  by  (P  '  P  ' )  and  (P  "  P  ").   Furthermore  for  each 

point  Px,   ,   >  on  the  three  dimensional  space  there  is  one  and  only  one 
V,a  ,  b ,  c  ) 

set  of  points  Px'  and  Px"  in  the  reference  plane  such  that: 

a.  Px1  lies  on  the  Pj'Po'  line 

b.  Px"  lies  on  the  P  "P  "  line 


plane; 


c.    If  g  is  the  intersection  of  OG  and  P1  Px"  on  the  reference 


|g  Px*|  =  b  ,   |g  Px"|  =  a  ,   |0g|  -  c  . 
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5.  Planes  on  the  Reference  Plane 

A  plane  is  represented  by  two  of  its  lines  or  three  of  its  points 
in  the  3-D  space.   So  it  needs  four  straight  lines  or  six  points  in  order 
to  be  defined  in  the  reference  plane.   Notice  the  exceptions  when  the 
actual  plane  is: 

a.  Normal  to  the  Or  axis.   In  this  case  all  the  four  images 
of  any  two  lines  selected  to  represent  it,  will  coincide  and  will  be  nor- 
mal to  the  OF  axis. 

b.  Represented  by  its  intercept  lines  on  the  OBr  and  OAT 
planes.   In  this  case  two  straight  line  images  will  suffice  and  the 
other  two  will  coincide  with  the  Of    axis. 

6 .  Lines  Other  Than  Straight,  on  the  Reference  Plane 

Much  of  the  power  of  the  descriptive  geometry  representation  is 
lost  when  one  deviates  from  straight  lines  or  planes.   Fcr  example,  the 
intersection  of  a  straight  line  and  a  plane  can  be  determined  by  a 
standard  procedure,  but  of  a  curved  line  and  a  plane  the  intersection 
can  be  determined  by  trial  and  error  only. 

However,  any  continuous  line  in  actual  space  will  give  contin- 
uous images  in  the  reference  plane  and  vice  versa. 

So  if  one  knows  that  the  images  (P  '  ,  P  "),  (P  '  ,  P"),  (P'',  P") 
and  (P, ',P ,")  as  shown  in  figure  5-3b  projected  by  points  P, ,  P„ ,  P^ ,  and 
P,  which  formed  a  continuous  line  in  space,  then  one  can  anticipate  that 
the  projections  of  the  actual  line  segment  would  more  or  less  look  as 
shown  in  figure  5-3b  and  the  approximation  of  the  actual  curve  shown  in 
figure  5-3a  would  be  more  accurate  if  more  points  at  closer  distances  were 
available. 
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B.    GENERAL  THREE  VARIABLE  PARAMETER  SYSTEM  WITH  TIME  DELAY- -A  CONCEPT 

1.    A  simple  extension  of  the  case  of  a  system  with  two  parameters, 
entails  that  for  the  case  of  three  variable  parameters  A  (Alpha),  B  (Beta), 
T  (Gamma),  the  general  characteristic  equation  would  look  like: 

ic^s)  +  e"TsC2(s)J  +  b[di(s)  +  e"TsD2(s)|  +l^E1(s)  +  e"TsE(s)]  + 
Ab[fi(s)  +  e"TsF2(s)]  +  ArUl^s)  +  e"TsH2(s)]  +  Blfl  (s)  +  e"Tsi2(s)]   + 
I^K^s)  +  e"TsK2(s)]  +[li(s)  +  e~TsL2(s)]   =   0  (V-l) 


AB 


2.    The  solution  for  the  T  and  UL'   curves  for  this  system  must  follow 

n 

the  following  steps: 

a.  First  derive  two  sets  of  equations  (separating  once  again 
reals  and  imaginaries)  in  terms  of  £  and  UU  ,  e.g. 

A^(tVn)  +  B61<g,u)n)  +  re1(g,uJn)  +  AB01(C,UJn)  +  AFt^CO  +  BTcp^C,^)  + 

AfiTn^C,^)  =   -^(C,^)     ,   and                             (V-2a) 

A^2(C,i'n)  •+  B62(C,U)n)   +  re2(C,Xn)   +  ABi2(C,Xn)   +  ATT]2<£,u>n)   +  BTcp^,^)   + 

ABT>2(S,U)n)  =     -§2^»V                                                                                            (V-2b) 

b.  Follow  the  same  steps  as  in  the  two  parameter  case  by: 

(1)  Giving  actual  values  to  UU   and  C    to  form  the  certain 

n 

"constant  performance  characteristic  curve,"  therefore  making  Ot   ,6  , . . . 
and  a    ,   6   ...  constant  for  each  "£,x  "  pair. 

(2)  Solving  the  two  equations  of  the  form  of  (V-2)  but 
with  constant  coefficients,  for  the  three  unknowns  A,  B,  F,  by  applying 
some  restriction  on  the  solution,  e.g.  by  watching  a  certain  plane  cut  of 
the  A,  B,  T  space. 

c.  Put  the  partial  plots  together  to  reconstruct  if  possible 
the  whole  picture  in  the  A,  B,  Y   space. 
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3.  This  idea  of  taking  cuts  on  the  three  parameter  space  is  not  new 
and  in  fact  solutions  by  quantizing  one  of  the  three  parameters  (which  is 
equivalent  to  subsequent  cuts  normal  to  one  of  the  three  parameter  axes), 
have  been  successfully  presented.   The  only  difference  is  that  this 
character  of  analysis  has  a  trial  and  error  characteristic  which  of  course 
is  reduced  if  the  physical  limitations  of  one  of  the  variables  are  more 

or  less  narrow. 

4.  However,  by  taking  enough  cuts  with  planes  passing  through  one 
of  the  coordinate  axes  of  the  parameter  space,  one  can  solve  an  equi- 
valent number  of  two  dimensional  problems  which  can  be  displayed  on  the 
same  reference  plane  by  use  of  simple  descriptive  geometry  ideas. 

5.  These  ideas  can  be  very  easily  grasped  for  the  case  of  a  sys- 
tem simpler  than  the  general. 


C.    A  SYSTEM  WITH  LINEAR  CHARACTERISTIC  EQUATION  IN  A,  B,  AND  T 

1.  Consider,  now,  equations  (V-2)  for  the  special  case  that  AB, 
AT,  BF,  ABF  do  not  exist.   In  this  case: 

a^ccv  +  B6i^'V  +  rei(£'V  =  _H(C^n) 

A*2(C^n)  +  B62(£,uM  +  re2(£,U)n)  =  -§2(C^n)  (V-2a) 

2.  Consider,  now,  that  one  wants  to  trace  the  "constant  uu   loci" 

'     '  n 

for  uu  =  uu   with  either  £   =  f,+£      or  f   =  £n  ,  F„,  .  ..,  f      for  a  contin- 
n    n^  s    sl  sn    s   sl  ^2  ^n 

uous  or  quantized  solution  (this  makes  no  difference).   That  is,  the 
following  set  of  equations  must  be  solved,  "n"  times  over  the  "£"  range 
desired: 


Aa2(c,.un)  +  bp^c,^)  +re1(c,u)n)  =  -52(C,uin) 


(V-2b) 


So  one  must  solve   "n"  times  for  the  intersections  of  two  planes  in  the 
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ABr  space.   In  other  words  the  constant  (JJ  =  UU    loci  is  a  set  of  "n" 
r  n    n\ 

straight  lines,  with  n  going  toCO  if  £  is  varied  continuously  in  the  range 

3.  But  by  letting  B  =  0,  for  any  arbitrary  £  =  £    from  the  above 

3  r* 

range,  in  both  of  equations  (V-2b)  and  solving  for: 

*"1«ar>'"m)+rel«ar.V>--?lCar'%l> 
to2«ar-V)+rel«ar>V)  =  -§2«ar'%1> 

one  gets  a  point  P,  in  the  OAF  plane  which  belongs  also  to  the  desired 

solutions  of  V-2b  for  C   =  f      ,  .       (see  figure  5-4a) . 

-    ^arbitrary 

Analogously  for  A  =  0,  one  can  get  a  point  P_ ,  on   the  OBr 
plane,  which  point  will  also  belong  to  the  desired  straight  line  solution 
of  (V-2b)  which  after  all  is  completely  defined. 

4.  So  this  solution  can  now  be  projected  in  the  descriptive  geo- 
metry reference  plane  (Figure  5-4b) ,  and  mapped  by  its  two  straight  line 
images.   These  two  images  define  any   point  that  happens  to  have  a  pair 

of  performance  characteristics  UU   =  1    ,  and  £  =  t      ,  ,  •   C   ,  £  (Ci"="C  )  L 

n    nl  arb  i_  arb   M  ^n  J 

for  example  a  point  Px  in  figure  5-4b  with  A  =  1,  B  =  1,  F  =  4  would  meet 
those  requirements.   Some  work  on  this  has  been  presented  by  Thaler  and 
Cadena . 

5.  The  fortunate  fact  to  note  though,  is  that  one  can  keep  track 

and  interpret  the  meaning  of  the  two  points  P, ,  P9  as  part  of  a  simple 

infinite  solution  for  UU   =  !x   ,    and  C   -  C      ,  ,  by  simply  locating  the  points 

n    nl     *   ^»arb   J  r  J 

P„ '    ,  and  P  ''   on  the  reference  plane,  since  P  ''   and  P.  '    can  be 

located  by  drawing  the  normals  to  Or  for  each  of  them.   One  should  of 

course  indicate  which  pair  of  'JJ  ,  and.  C      ,  was  used  for  this  set  of 

v  nl     sarb 

points . 

6.  But  now  looking  back  again  in  the  whole  process  one  recognizes 
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that  what  is  necessary  to  solve  this   problem  is  to  solve  simultaneously 

for  each  set  of  values  of  uu  =  UU  .  and  F  =  F   ,  _     ,  a  two  dimensional, 

n    nl     *   ^arbitrary 

linear  parameter  plane  problem  in  A,F  and  B,r  and  to  locate  also  simul- 
taneously the  two  corresponding  points  on  the  same  planar  plot.   This 
process  will  continue  regularly  in  order  to  complete  two  individual  planar 
grids  of  constant  uu  and  constant  £  curves  as  shown  in  Figure  5-5. 

7.  It  is  important  to  plot  each  of  the  two  grids  with  different 
locus  indicators,  eg   for  the  OBF  and    for  the  OAF  projections. 

8.  If,  after  the  two  parameter  plane  grids  are  plotted,  one  wants 

information  about  the  uu  =  uu  ,  ,  F  =  F,  ,  A,B,F  space  locus,  the  procedure 

n    nl      M 

is: 

a.  First  have  in  mind  that  the  plot  obtained  lies  in  the 
reference  plane. 

b.  Second,  locate  the  two  points  for  each  grid,  corresponding 

to  the  uu  =  uu  ,  ,  F  =  C-,    loci, 
n    nl   J    ^1 

c.  Draw  the  two  normals,  on  or  to  locate  the  associated 
second  projection  for  each  image. 

d.  Draw  the  lines  il    and  iu   which  are  the  images  of  the  de- 
sired locus. 

e.  In  order  to  read  one  of  the  infinite  sets  of  A,B,r  values 

that  give  UU  =  UU  .  and  F  =  F,  : 
n    nl     *   M 

(1)  First,  select  a  F  value,  say  F  =  Og 

(2)  To  get  the  other  two  parameter's  values,  sketch  a 
line  normal  to  gOF,  at  the  point  g.   Then:   A  =  ga,  B  =  gb,  the  intercepts 
of  the  line  bga  on  the  projections  X"  and  X1,  respectively. 

9.  An  example  of  an  actual  simple  linear-characteristic  equation 
system  follows. 
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D.    EXAMPLE  PROBLEM 

1.    The  closed  loop  transfer  function  for  the  system  in  figure  5-6 
is: 

c r'(s-i) _ 


R    (s2+  3s  +  2)  +  r'(s-l)(A*s  +  B*) 


(s  -1) 


T(s2  +  3s  +  2)  +  (s-l)(A's  +  B  +  1) 


where 


r  =  jrr      ,  and   B  =  B'  -  1. 


r 

2.  So  the  characteristic  equation  for  this  system  being  of  the  form 
previously  discussed  is: 

A'(s2  -  s)  +  B(s  -  1)  +  T(s2  +  3s  +  2)  +  (s  -  1)  =  0     (V-3) 

3 .  Practical  Considerations  on  Applying  Parameter  Plane  Methods 
At  this  stage  it  is  of  importance  to  notice  that  a  digital 

computer  program  for  solving  and  drawing  the  constant  UU   and  constant  £ 
loci  is  absolutely  necessary.   One  could  use  either  a  program  based  on 
the  equations  of  paragraph  E  of  this  section  (V) ,  or  could  use  twice  a 
program  designed  to  solve  and  draw  for  the  case  of  two  variable  para- 
meters, and  then  superimpose  the  two  drawings  (one  for  an  AOF-plane  cut 
and  the  other  for  a  BOT-plane  cut).   But  even  if  one  uses  the  two  vari- 
able parameter  plane  program,  difficulties  may  arise  in  selecting  the 
proper  range  of  A,B,r  (that  is  the  scale  of  the  graph)  basically  for  the 
desired  range  of  "d"  and  "£" .   However,  even  if  one  or  two  curves  are 

drawn,  the  task  then  is  even  simpler. 

*2 

To  get  a  first  rough  estimate  of  the  appropriate  range   (under 


*2 

Refer  to  the  conclusions  (section  V-G)  for  comments  regarding 

root  loci  methods  versus  parameter  space  approach,  and  also   note  that 
the  solutions  by  the  root  loci  techniques  in  subparagraphs  "a"  and  "b" 
to  follow,  are  just  particular  arbitrary  solutions  and  are  not  likely  to 
be  the  desired  solution. 
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the  assumption  that  some  test  runs  with  the  computer  program  of  Appendix 
D  will  follow),  one  can  proceed  in  many  ways.  Three  of  them  are  as  fol- 
lows: 

a.  The  Case  of  Forward  Loop  or  Feed-Back  Compensation  with  the 
Forward  Gain  as  One  of  the  Variable  Parameters 

(1)  Draw  the  open  loop  root  locus  with  rio  compensation 
unit  included 

(2)  Acknowledging  the  fact  that  except  for  the  forward 
gain,  the  rest  of  the  variable  parameters,  once  fixed,  define  possible 
compensation  poles  and  zeros,  try  to  find  a  combination  of  these  para- 
meters that  would  change  the  open  loop  root  locus  to  be  stable  for  the 
desired  working  frequency  with  a  £  approximately  equal  to  0.707. 

(3)  As  a  very  quick  check  on  the  approximate  forward 
gain  range  apply  the  Routh-Hurwitz  stability  criterion  for  the  closed 
loop  characteristic  equation.   Sometimes  this  is  not  enough  and  one.  should 
draw  the  root  locus. 

(4)  The  rough  scale  to  be  used  for  the  first  time  should 

*3 

be  as  a  rule  of  thumb    ten  to  a  hundred  times  the  values  of  the  com- 
pensation variables  selected  and  the  maximum  stable  forward  gain. 

b.  Compensation  (Forward  Loop  or  Feed-Back)  with  a  Specified 
Forward  Gain 

(1)  Insert  a  variable  gain  parameter  in  place  of  the 
specified  forward  gain  value. 

(2)  Repeat  the  procedure  for  the  previous  case,  but  ac- 
cepting the  resulting  scale  if  the  specified  gain  is  as  a  rule  of  thumb 


*3  df 

The  rule  of  thumb  proposed,  might  of  course  not  work  for  -rr 

very  large   or  very  small,  and  therefore  some  change  should  then  be 

necessary  depending  of  the  particular  system. 
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twenty  percent  less  than  the  maximum  stable  gain  obtained.   Otherwise, 
try  a  new  combination  of  the  compensation  variables.   The  reason  of  this 
procedure  is  to  assure  an  A,B  pair  which  would  give  under  the  operating 
"k" ,  at  least  one  complex  pole  pair  in  the  region  of  £  =  (.3  t  .8),  be- 
cause the  possibility  of  not  getting  a  constant  £  curve  is  reduced  when 
this  pair  is  forced  in  the  corresponding  s-domain  area,  for  "£"  around 
0.707. 

c.    The  Closed  Loop  Characteristic  Equation  Technique 
'     When  only  the  closed  loop  characteristic  equation  is 
given,  and  one  has  no  idea  about  the  actual  system,  the  following  pro- 
cedure could  be  followed: 

(1)  Obtain  three  Root-Locus  Plots,  each  for  one  of  the 
three  parameters,  as  the  only  variable  and  the  remaining  two  set  equal 
to  a  convenient  constant  (say  unity). 

(2)  For  a  certain  range  of  "£"  and  "Jj  "  desired,  define 
the  corresponding  range  of  each  of  the  variables  in  each  of  the  root- 
locus  plats.   This  will  give  an  indication  of  the  scaling  of  the  problem 
to  start  with. 

One  should  also  have  in  mind,  once  engaged  with  a  time 
delay  problem,  to  set  the  time  decay  constant  equal  to  zero,  to  apply  the 
method  discussed  above. 

4.    Practical  Considerations  on  Applying  Parameter  Methods  of 
Paragraph  C 

Considering  equation  (V-3),  and  having  in  mind  the  discussion  of 
paragraph  C,  one  comes  to  the  conclusion  that  if  a  BOF  cut  of  the  para- 
meter space  is  tried  for,  the  particular  system  of  figure  5-6  (i.e.  A=0) 
is  physically  restricted  from  having  any  £,  since  the  roots  of  the 
GH  function  will  all  lie  on  the  a  axis  in  an  "s-plane"  representation. 
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So  there  is  a  physical  limitation  in  applying  the  technique, 
which   one  can  overcome  by  a  change  of  variable  A1  =  A  +  1,  thus  main- 
taining the  tachometer  feedback  of  the  system  for  the  case  of  A  =  0.  In 
this  case,  equation  (V-3)  becomes: 

A(s2  -  s)  +B(s  -1)  +  T(s2  +  3s  +  2)  +  (s2  -1)  =  0   (V-4) 

5.    Having  in  mind  all  the  practical  considerations  above,  one  can 
proceed  to  specify  an  approximate  set  of  A,  B,  T  range  for  a  first  try 
of  the  parameter  space  solution. 

a.  From  figure  5-7  (following  the  procedure  section  V-D3) , 
one  can  see  a  good  value  for  the  compensation  zero  to  be  at  a  =  2 , (or 

0"  =  3 ,  etc.)-   This  can  result  for  A'  =2,  B'  =  -4  and  equivalently  for 
A  =  1,  B  =  -5. 

b.  Since  the  Routh  Hurwitz  criterion  is  going  to  give  an 

expression  of  F'    <  M   .  ,  for  M   .   >  0,  it  follows  that: 
max    gam       gain 

r       =  ( — - — )  >  m    . 

max    p,        gam 
max 

for  stability,  and  this  is  not  enough.  (Root-locus  is  needed). 

c.  For  A  =1,  B  =  -5,  the  second  order  system  of  equation 
(V-4) ,  becomes : 

s2  +  3£-2}s 
s    +     r+2 

It  follows  then  that  for  a  C    =  1//2-  and  JU  =  /l~        a  T  =  10  is  needed 

*  n 

d.  So  as   a  first  try  the  following  ranges  should  be  tried 
out: 

(i)  |a|  s  io 

(2)  |B|  ^50 

(3)  |r|   ^  20 
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6.    In  figure  5-8,  the  results  of  the  parameter  space  analysis  are 

shown.   In  order  to  show  that  the  results  are  valid,  select  the  parameter 

space  loci  for  (a1  =  1.0,  £  =  0.7),  represented  by  its  two  images  (the 

dash-dotted  lines).   An  arbitrary  set  of  values  for  this  loci  is, 

(A  =  8.375,  B  =  2.32,  V   =  12.0), 

which  gives  a  second  order  polynomial,  (equation  V-4) : 

2   29.625    20.680 
S  +  21.375s  +  21.375 

with:  uu  =  .985 

,      n 

C   =  .706 

For  a  second  demonstration  consider  the  loci  for  (uu  =  1.6,  F  =  .6).   An 

n 

arbitrary  set  of  values  for  A,  B,  T   on  this  loci  is, 

(A  =  -.718,  B  =  -3.5,  r  =  3.0), 

giving  a  second  order  polynomial , (equation  V-4): 

2   j^.218    8.500 
s  +  3.282s  +  3.282 

with:  uu  =  1.62 

n 

£   =  .603 
These  results  turn  out  to  be  fairly  accurate  and  the  theory  presented 
works  in  practice. 

E.    SOLUTION  TO  THE  GENEPvAL  PROBLEM 

1.  Consider  again  equation  (V-2) .   The  discussion  in  section 
(V-B)  indicated  the  necessity  for  an  infinite  number  of  cuts  to  be  done 
for  the  complete  parameter  plane  solution  in  three  dimensions.   However, 
depending  on  the  desired  degree  of  accuracy  a  finite  number  of  cuts  will 
suffice  for  practical  purposes. 

2.  The  solution  for  a  case  with  more  than  two  cuts  will  be  demon- 
strated in  section  5. 
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3.  The  equations  for  the  general  cut  through  the  T-axis  will  be 
derived  by  letting: 

A  =  XB  (V-5) 

as  the  plane  of  the  cut  implies. 

4.  By  substituting  equation  (V-5)  into  equation  (V-l),  the  fol- 
lowing set  of  two  equations  results: 

B(A^1  +  61)  +  Tel  +   B2\0L  +  Br^  +  cp  )  +  B2r\M,   =  -?       (V-5a) 

b(\<*2  +  &2)  +  re2  +  b2\02  +  Br(\Ti2  +  cp2)  +  b2hv2  =  -5  (v-5b) 

Therefore: 


r  = 


-§1    -Q,u1  +  61)B   -^B' 
el  +B(\ti1  +  cpx)   +B  X|J.1 


(V-6a) 


and 


BXA^   +  B     i\     (A.+A,)+\(Ar+A      )         + 

2  r  2  n 

B     I  \   A  +  \  ( A        +  A,+AJ+A_)+AS  + 

L       rp  cpa         n&         e0         (j,s  cpo     J 


Tpr  *  cpa  r|6 

B     J  \(A       +  A  _)   +  (A  _  +  A  .  )~|  +  A 

L       ea  tic,  cp^         e6   J          ( 

where  in  general 


e§ 


=  0 


(V-6b) 


A 


xy 


Xl   yl 


x2  y2 


Xiy2  "  yiX2 


5.    So  in  order  to  define  a  certain  desired  "(JU  ,  f   "  point  in  the 

t'  st   F 

A,B,F  space,  one  should  follow  steps  as  follows: 

a.  Define  the  coefficients  Oi    ,  (3  ,  ...,  §,  and  »■,  £_,...§„, 

which  become  constants  for  UJ  =  UU  .   T   -  C    •   It  is  assumed  of  course  that 

t  *     t 

the  time  delay  in  the  system  is  fixed. 

b.  Solve  for  B  from  the  equation  (V-6b).   A  digital  computer 
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routine  will  be  necessary  for  this  purpose.   However,  the  conventional 
solution  of  a  quartic  equation  is  listed  as  Appendix  E,  for  easier  ref- 
erence.  Care  should  be  taken  to  exclude  any  complex  roots. 

c.  Then  T   will  be  given  by  using  (V-6a),  for  each  admissible 
B-root. 

d.  If  track  of  the  cut-plane  is  kept  the  other  projection  is 
not  necessary,  since  A  =  \B  (V-5),  and  therefore  for  each  point  in  the 
BT  grid  that  will  be  generated  by  the  procedure  of  subparagraphs  5a 

5c,  there  will  correspond  another  one  as  shown  in  figure  5-9. 

In  figure  5-10a  the  grids  of  constant  £  and  uu  curves 

obtained  for  three  different  cuts  are  shown.   Notice  that  one  projection 

for  each  cut  is  sufficient. 

To  construct  an  approximate  constant  JU  ,  £      curve  one 

n  ^n 

must  select  the  points  P' ,  P '  .  P~   on  the  corresponding   grids  (  one 
from  each  available  grid),  and  by  joining  them  by  straight  line  segments 
obtains  one  of  the  two  required  approximate  projections. 

In  figure  10-b  this  projection  is  shown,  with  the  con- 
struction grid  eliminated.   The  second  projection   of  the  constant  UU  , 
£   curve  can  be  approximately  defined,  by  locating  the  projections 
P",  P",  P"  from  equation   (V-5),  with  X    =  X    ,  X    ,  X    ,  etc. 

F.    EXAMPLE  ON  A  SIMPLE  DESIGN  PROBLEM 

Consider  the  hypothetical  system  of  figure  5-11,  and  suppose  that 
the  operational  requirements  of  the  system  are  met  if  the  operating  point 
is  on  =  .4  rad/sec  and  f      =  .5. 

Consider  also  that  the  cost  for  building  this  system  was  found  to 

-Ts 
be  C  =  M»r  where  M  is  a  big  number  and  T  is  so  small   that  e    =1 

(in  order  to  simplify  the  algebra).   There  exists  no  information  about 
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BLOCK  DIAGRAM  FDR  A  SYSTEM 
WITH  TIME-  PEL  A  Y 
(PRODUCTS  INCLUDED) 
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the  operating  cost  of  the  system  as  a  function  of  A,  B,  F,  so  this  aspect 
will  not  be  considered,  but  of  course  the  stability  is  presupposed. 

It  is  desired  to  select  A,  B,  F  that  will  simultaneously  give  UU  , 
£   as  specified  and  minimum  possible  cost. 

It  is  noted  that  A  physically  represents  the  positive  gain  of  *an 
amplifier,  therefore  there  exists  the  restriction  A  >  0. 

SOLUTION 

1.    The  closed  loop  gain  of  the  system  turns  out  to  be: 

_C  =  r(s+3)(As  +   l)e"Ts 

R       [s4+    (2.828  +B)s3+   (4+2.828B)s2+  4Bs]  +e"Ts[Ar(s2+  3s)    +F(s+3)1 

So   the  characteristic   equation   for   the   general   cut-plane  A  =  XB,    and 
T  -•  0  will   become: 


T4  3  2~>  T    3  2  ~\ 

CE(s)   =  Is     +  2.828s     +  4s  J  +  b|  s     +  2.828s     +  4s  !   +  F 


s+3   + 


r  2     "i 

BFi  X  s  +  3\  s   =  0 . 

L        J 

2.  At  this  point  one  must  make  a  decision  on  selecting  the  number 
of  cuts  that  will  suffice  the  specific  application.   For  this  case  con- 
sider eight  cuts,  corresponding  to: 

\   =  o,  X   =  ±%,  X   =   ±1,  X   =  ±2,  X   =  oo 

In  table  5-1  the  resulting  CE(s)  are  tabulated. 

3.  By  applying  the  method  of  section  VE,  the  curve  for  constant 

(uu  =  .4  and  £  =  .5),  was  obtained  and  displayed  in  figures  5-12a,  5-12b, 
These  figures  have  been  plotted  according  to  the  data  of  table  5-2. 

In  figure  5-12a  the  interesting  situation  for  the  stable  com- 
pensation of  the  system  under  design  is  displayed. 

4.  According  to  the  available  information: 

a.    Gain  Upper  Limit:  F=0.217  units  (stable  operation  on 


specified  UU  ,  f) 
n.  s 
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b.  Gain  Lower  Limit:   T  >  0.0.   Nevertheless  it  was  possible 
to.  get  a  minimum  gain  of  .184  units  for  the  lowest  B  parameter  employed. 

c.  Compensation  Pole: 

0.00284  <  B  <  0.47  in  the  range  T  =  (0.184  t  0.217) 

d.  Tachometer  Feed-Back: 

0.000467  <  A  <  2.84  in  the  same  T  range.   So  it  follows 
that  the  system  is  more  sensitive  in  "B"  than  in  "A"  variations  (for  this 
specific  T  range). 

5.  So  in  order  to  satisfy  the  requirement  of  as  small  gain  as 
possible,  imposed  by  the  specifications  and  having  the  information  of 
figure  5-12a,  one  should  choose  of  course  the  set: 

(A,  B,  O   =  (2.84,  .00284,  .184). 

*4 

But  if  it  was  possible  to  pay  9  percent  more   the  solution 

(A,  B,  T)2  =  (.115,  .25,  .20) 
is  much  more,  advantageous  and  takes  care  of  possible  component  drift. 

6.  By  substituting  the  set  (A,  B,  T)    =  (.115,  .25,  F)  in  CE(s), 
(for  T=0),  one  gets: 

CE(s)  =  '  s4  +  3.078s3  +  4.707  s2  +  1.00s!  +  H"  1.15s  +  4.445s+3J  -  0. 

the  root-locus  plot  of  which  gives  (ju  X)  .   .    =  (.3986,  .477)  for 

n  -  dominant 

r  =  .2,  which  is  very  close  to  the  desired  specifications   (see  figure 
5-12c). 

7.  It  is  of  interest  to  note  what  the  constant  a  =  0  and  £  =  0 
surfaces  look  like. 

a.    From  CE(0")  for  s  =  (a  =  0)   the  a  =  0  surface  is  the 
T  =  0  plane. 


r  -r 

-/  /\t>        /        1 

Since  the  cost  function  is  C  =  MT  it  follows  that  — -= -p -.087 

Cl   ri 
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b.  From  CE(jJu)  for  the  s  =  jju  surface  one  gets  the  following 
set  of  equations. 

u)4  -    (4  +  2.828B)dJ2   -  Aruu2  +  3T  =  0 

-(2.828  +  B)i'3  +  4BUJ  +  3ATuu  =  0 

c.  By  simultaneously  solving  the  O"  =  0  and  £  =  0  equations 
the  equation  for  the  cut  of  those  two  surfaces  is: 

B2  +  2.828B  +4=0 
which  gives: 

b  =  -  /i  ±  j  ji 

Therefore  there  exists  no  intersection  of  the  a  =  0  ,  £  =  0  surfaces. 

G.    CONCLUSIONS 

From  the  theory  presented  in  this  section   theoretically  and  shown 
to  actually  work  for  the  simple  system  of  figure  5-11,  the  following  con- 
clusions can  be  obtained: 

1.  It  is  possible  to  construct  the  constant  "J.'-£"  curves  for  a 
three  variable  parameter  system.   Care  must  be  taken  in  order  not.  to 
exceed  the  physical  limitations  of  the  system,  e.g.  taking  cuts  that 
would  perhaps  destroy  the  feed-back  ot»  forward  gain  of  the  system. 

2.  In  the  case  where  the  variable  parameters  show  in  a  nonlinear 
fashion  in  the  system,  the  constant  "uj-£"  curves  have  no  straight  line 
projections . 

3.  Care  should  be  taken  to  differentiate  between  different  seg- 
ments of  the  solution.   It  is  of  importance  to  define  which  of  them 
provide  stable  and  which  unstable  operation.   Until  a  well  grounded 
stability  theorem  is  developed  in  the  three  dimensional  space,  this  task 
which  was  easy  for  the  system  of  section (V-Fjmight  prove  impossible 
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for  a  complicated  situation.   In  this  case  the  method  of  this  section 
will  be  handicapped. 

4.  Another  thing  to  notice,  is  that  for  the  nonlinear  case,  (A,B,r 
products),  it  is  possible  to  come  up  with  a  curve  which  is  more  or  less 
normal  to  the  T-axis.   If  this  happens  it  is  not  easy  to  read  the  appro- 
priate A,B,r  sets  in  one  hand,  but  sensitivity  on  the  T   parameter  is 
proved.   So  if  one  wants  to  operate  under  the  conditions  imposed  by 

the  curve,  one  must  set  F  =  F      and  then  possibly  one  can  reduce  the  prob- 
lem to  a  two  dimensional  one. 

5.  Sensitivity  analysis  of  a  system  is  possible  for  a  multivariable 
cor>trol--system,  by  applying  methods  described  in  this  section.   This  is 
in  contrast  to  the  approach  of  using  a  family  of  root  loci  where  the 
method  of  sensitivity  analysis  is  not  so  direct. 

6.  Another  great  advantage  of  the  parameter  space  approach  over  the 
root  loci  method  is  that  system  optimization  through  adjustment  of  para- 
meters is  more  direct,  and  illustrative. 

7.  Of  course  in  order  to  arrive  at  the  desired  result  a  vast 
amount  of  analysis  was  necessary,  and  it  is  a  disadvantage  which  of 
course  is  removed  when  a  computer  program  similar  to  the  one  in  Appendix 
D  is  available.   So  indeed,  the  method  is  "digital-computer"  oriented. 


123 


VI.   AN  ACTUAL  ENGINEERING  APPLICATION 

A.  GENERAL 

1.  In  this  section  the  roll  angle  control  system  for  a  vertical 
take-off  aircraft  is  studied. 

2.  In  a  post-design  analysis  of  this  system,  (which  was  initially 
designed  by  the  root-locus  methods  and  thereafter  operated  successfully), 
attempts  to  locate  the  actual  operating  point  inside  the  two  variable 
parameter  stability  region  were  originally  unsuccessful.   This  implied 
that  the  stability  criterion  of  chapter  IV  was  invalid,  and  since  the 
derivations  of  the  stability  criterion  was  rigorous,  further  study  of 

the  particular  case  is  needed. 

3.  In  the  subsections  to  follow,  the  system  is  first  considered 
from  an  elementary  aspect  and  since  it  needs  compensation,  its  possible 
modes  are  defined. 

Then  a  parameter  plane  design  follows  and  having  in  mind  the 
results  of  the  stability  analysis t   an  appropriate  operating  point  is  se- 
lected.  Time  domain   simulations  follow  at  last  to  insure  the  correctness 
of  the  parameter  plane  results. 

B.  THE  SYSTEM 

1 .    Brief  Physical  Description 

a.    The  system  is  presented  in  block  diagram  form  in  figure 
6-1. 

In  the  component  analysis  of  this  figure  one  can  observe 
seven  variable  parameters.   But  since  the  servo  constants  have  been  pre- 
set to  their  values  of  table  6-1,  the  only  variable  parameters  remaining 

are  K  ,  K  ,  and  T  (which  is  bounded  within  limits). 
P   0 

124 
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TABLE  6-1 
VALUES  OF  CONSTANTS  USED  IN  THE  SYSTEM  OF  FIGURE  6-1 

CONSTANT   VALUE 

T  .15  t  .30 

K  1/8.1  =  .1235 

D  .488  x  10"3 

E  .4719  x  10"1 

F  2.329 
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b.    Since  the  permissible  range  of  variation  on  T  is  small,  it 
is  reasonable  to  start  with  the  minimum  expected  value  of  T,  (0.15),  since 
the  results  of  section  IV  show  that  this  should  give  the  maximum  stability 
region  in  the  A,B  plane. 

2 .  Objectives 

a.  To  actually  define  appropriate  K  ,  K.  values  for  stable 
operation.   In  doing  so  one  should  have  in  mind  that  zetas  around  0.7 

and  frequencies  from  1  to  3  rad/sec  are  desirable   from  a  "pilot-perform- 

*5 

ance"  point  of  view. 

b.  To  check  and  find  if  the  operating  values  of  K  ,  K.  lie 
in  the  parameter  plane  stability  region  in  order  to  check  on  the  dis- 
crepancy reported  (section  VI  A).    These  values  are: 

(K  ,  K  )  =  (95,  190) 

P    0 

3 .  Systems  Transfer  Characteristic 


K 

&,   N    (.305)  D  e   °Ls"  +  .444s  +  .488.J  ,__.  ,  N 

*(s)  =  -> ^ (VI-1) 


"TSI  2  .488.J 


where: 


and 
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A'P  (s)  =  A  s3  +  .444s2  +  .488s  +  O.oj 

r-    3     2  "1 

B'PT3(s)  =  B  •  0«s   +  s   +  .444s  +  .488  I 

B         L  J 


A  =  (-325)K  K  =  77.2  K    =»   Tachometer 
Dp         p 

B  =  (»305)K    _    2      ^   Proportional  Control 

D     0         0         v 

F(s)  =  s6  +  98.4696s5  +  4941.1680s4  +  8978.082s3  + 
5086.9961s2  +  3162.6011s  -  147.59 


■'5 


These  were  established  by  the  original  system  designer. 
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4.    The  Systems  Characteristic  Equation  is: 
CE(s)  - 


'  s6  +  98.4696s5  +  4941.1680s4  +  8978.082s3  + 


5086.9961s2  +  3162.6011s  -  147. 59j  + 


e"rs  A  Is3  +  .444s2  +  .488s  +  O.o]  + 
e~Ts  B  [  0's3  +  s2  +  .444s  +  .488 
=  0  (VI-3) 


C.    POSSIBLE  COMPENSATION  MODES 

1.  By  looking  closer  to  the  feed-back  component  of  the  figure  6-1 

system,  we  can  see  in  general  that  its  effect  on  the  open  loop  transfer 

function  is  the  insertion  of  a  compensation  zero  "z  ",  and  by  controlling 

the  values  of  the  variable  parameters  K  and  K_,  one  can  cause  various 

P       0 

compensation  modes. 

2.  Noting  also  the  results  of  section  II,  where  it  was  observed 
that  though  the  shape  of  a  root- locus  was  symmetrically  deformed  towards 
the  jJJ-axis,  its  basic  form  remained  unchanged  with  varying  delay,  one 
can  search  for  the  modes  mentioned  above  using  the  delay-free  version  of 
figure  6-1  without  loss  of  generality. 

3.  In  figure  6-3  the  open-loop  root  locations  of  the  system  are 
shown,  with  compensation  z   not  included. 

So  with  no  compensation  the  pole  P,  will  produce  instability  in 

the  desired   operating  range. 

Note  also:   first,  that  poles  Pc  and  P   located  at  a  natural 

j      6 

frequency  contour  of  69  rads/sec  with  a  £  =  .7,  are  too  far  away  to 

be  considered  at  the  desired  frequency  range  and  second,  that  zeroes 

z.  ,  z   restrict  root-travel  branches  leaving  poles  P„,  P.,  to  remain 

in  a  very  limited  range  in  the  s-plane,  so  that  they  force  the  continuous 

presence  of  a  complex  root  pair  with  a  set  of  (£ ,  i  )  =  (.43,  .7)  on  the 

average.   Theoretical  cancellation  would  of  course  happen  for  infinite 
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gain,  but  since  these  four  roots  are  very  close  together  in  pole-zero 
pairs,  the  influence  of  (P„,P. )  in  the  systems  transient  characteristic 
will  be  much  reduced. 

4.  Now  the  real  zero  "z  "  introduced  by  the  feed-back,  can  general- 
ly take  three  possible  relative  locations  with  respect  to  the  P, ,  P„  pair. 
as  follows: 

a.  z   to  the  Left  of  P_  . 

c  2 

This  will  cause  a  closed  form  of  the  root  locus  and  for  a 
proper  gain  will  provide  a  pair  of  complex  roots  in  the  L.H.P.  solely. 

In  this  case,  both  K  and  K,  must  be  positive. 

P       0 

b.  z   to  the  Right  of  P  . 

In  this  case  one  of  the  two  feed-back  gains  must  be  nega- 
tive. So  if  in  the  feed-back  expression 

K  s  +  K  ,  k  is  factored  out  of  the  parentheses 

p      0    p 


one  will  have: 


K 

■  K 


K  (s  +  -^) 


P 

which  implies  a  negative  gain  locus,  for  K  <  0. 

c.    z   Between  P_  and  P.  . 
c  2      1 

This  will  result  to  a  real  root  pair,  one  of  which  will 

be  in  the  right  half  plane,  for  "low-enough"  gain. 

5.    After  this  discussion, it  is  obvious  that  two  are  the  possible 

compensation  modes,  one  with  the  "z  "  in  the  L.H.P.  and  one  with  "z  "  in 

c  c 

the  R.H.P.   Note  that  the  operating  point  of  the  system  was  regulated  for 
a  Left  Half  Plane  compensation  with: 
(K  ,  K.)  =  (95,  190) 

P    0 

giving  a  zero  at  z   =  -2. 


130 


ROOT  LOCUS  FOR 
(Kp,Kf)=(95J30) 


LEGEND 


O,  x 


:    ROO T LOCUS 

:    ROOT,    LOG  AT 7 ON S 

FOR   &=  0 

:  P±  ,P2    F0£  G  =  7520 


Uu> 


i2 


s 


/ 


_  _  _  _<$_ _,._.*, Q.- 


3=-*i2    v  .  to 

t     G=2240 

-x  -««►--» — «*- x >•  6 

1/ 


^ 


\ 


^ 


♦ 


/ 


I  c  2 


Scale 

FOR   BOTH  dj'u) 
£  UNIT  P£R   INCH. 


Ul 


FIGURE  6-4 


ex 


131 


6.    In  figure  6-4a,  the  root  locus  of  the  system  is  plotted  for 

this  selection  of  (K  ,  K, ) .   This  was  done  as  usual  by  equating  the  open- 

P   0 

loop  transfer  function  to  minus  one,  and  thus  getting: 

[s  +  98.5s5  +  5056.3s4  +  8966.60s3  +  6164.65s2  + 
4905.43s  -  225.3  j  +  G  is3  +  2.443s2  +  1.371s  +  .973j  =  0 

where  G  was  the  variable  gain  parameter: 


,G  = 


r (.305)  x  K  x  K 
D 


"(Kp=95) 


=  59400  K, 


and  since  the  working  K  of  the  system  was  selected  to  be  .1235,  one  gets 
the  operating  point  to  be  as  follows: 


7320 


K 


.1235 


K 
P 

95 


0 

190 


B 


c 


7340    14680     .812 


uu 

j  /  n 
rad/sec 


1.65 


7.    To  illustrate  by  the  root-locus  technique  that  compensation  is 

possible  for  z   in  the  Right-Half  Plane,  one  can  select  (K  ,  K,)=(-95,95) 

c  P    0 

and  get  from  the  open  loop  transfer  function: 

Is  +  98.5s  +  5056.3s  +  8966.60s  +  6164.65s  + 
4905.43s  -  225.3  j  +  G*[  s3  -  .557s2  +  .044s  =  .487]   =  0 


where: 


(.305)  x  K  x  K  _ 


G'  = 


L 


D 


=  -59400  K 


(Kp  -  -95) 


So  for  positive  K,  as  in  this  cass  a  negative  gain  locus  is  resulting 

and  is  displayed  in  figure  6-4b.   In  this  case,  K     for  stability,  is 

r  J  °  max 

equal  to  .0965  (equivalent  for  G1  =  -5722),  which  is  a  restriction  for 
the  designer.   Notice  also  that  for  "£"  in  the  area  of  .707  or  greater, 
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the  natural  frequency  will  be  less  than  .55rad/sec. 

8.    With  this  information  about  the  system,  its  preliminary  investi- 
gation is  concluded  and  the  parameter  plane  design  follows. 

D.    DESIGN  IN  THE  PARAMETER  PLANE 

1.  By  using  the  computer  program  of  Appendix  D,  and  the  selected 
sets  of  constant  "£"  and  constant  "uu  "  curves  shown  in  table  6-2,  three 
sets  of  parameter  plane  grids  were  plotted  for  values  of  the  time  delay: 

(T1'  V  V  =  ('15'  ,2'  -3) 
These  results  have  been  plotted  in  figures  6-5a,  6-5b,  and  6-5c.   Quasi- 
linear  versions  of  the  maximum  root  regions  found  in  chapter  IV  have  been 
superimposed  for  easier  reference. 

2.  Having  in  mind  the  general  range  of  (V , '    )  desired,  (section 
VIB),  and  the  particular  (£  ,U)  )  =  (.812,  1.65)  that  resulted  to  be  the 
ones  preferably  selected  by  the  designers  of  the  system,  one  can  plot 
the  variation  of  the  (uu  =1.6  and  £  =  .8),  constant  (UJ  ,£)  curve  for 
varying  time. delay.   This  plot  and  its  expanded  scale  version  both  re- 
sulting from  data  of  table  6-3  are  shown  in  figures  6-6  and  6-6a.   Notice 
the  approximately   linear  relation  in  these  plots,  especially  within  the  . 
operating  time-delay  constants  T  €  [.15,  .3]  where  AT  and  displacement 

in  the  A,B  plane  remain  in  a  constant  ratio. 

TABLE  6-2 

Constant   "C"  r      =    .6     Co   =    -^     Co   =    •?     £/    =    .8     t,   =    .85 

Curve    Values  1  2  ^3  ^4  -5 

Constant  "uu"  „        ,OQ      .  _       n  „      ..  - 

n  ..   ,       uu.  =  .3  uu0  =  .688  uuo  =1.0  uu.  =1.3  uuc  =1.6 

Curve  Values     12         3        4        5 

This  fact  can  help  in  predicting  appropriate  A,B  values  for  opera- 
tion with  a  "time-delay  constant"  different  than  (T.,^,  T„)  =  (.15, .2, .3) 
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However  since  there  exists  no  limitation  on  the  lower  value 

of  the  time-delay  to  be  used,  or  in  other  words  since  T  .   =  .15  is  an 

nun 

admissible  delay  constant,  the  best  choice  of  A,B  values,  in  the  absence 
of  any  other  criterion  to  be  met,  is  the  set  for  T  =  .15: 
(A,B)  -  (5574,  9453) 

3.  After  having  the  first  complex  pole  pair  fixed,  one  should 
think  about  the  other  two  pairs. 

a.  In  the  case  the  system  is  operating  within  the  stable 
region,  one  should  take  into  account  that  the  complex  pair  at  (JO   =  69 
rad/sec  contributes  a  very  small  percentage  to  the  total  system  transient 
response  and  therefore  can  be  neglected. 

b.  If  the  second  pair  remains  near  the  (UJ  ,  T)  =  (-7,  -35) 

zero  pair,  in  the  case  now  that  finite  delay  is  present,  there  will  again 

occur  small  interference  on   the  transient  response  from  P~  ,  P..   But  in 

order  for  this  to  happen  one  should  be  able  to  observe  constant  (JJ  ) 

n 

curves  of  about  (.7)  all  over  the  stable  region  for  time-delay  constant 
in  the  range  from  (.15  v  .3).   To  show  that  this  is  true  the  parameter 
plane  curves  were  solved  for  an  average  time  delay  value,  (T  =  .225),  and 
plotted  in  figure  6-7.   So  by  observing  this  figure  one  can  see  the  ex- 
istence of  the  x   =  0.7  curve  spreading  out  across  the  A,B  plane  and 
therefore  present  for  every  (A,B)  selection. 

4.  As  far  as  the  solution  which  was  obtained  by  conventional 
means  (root-locus,  etc.)  one  can  observe  in  figure  6-6a  that,  for  T  =  0.0, 
the  set 

(A,B)  =  (7340,  14680) 
corresponding  to 

(K  ,  K  )  =  (95,  190), 
lies  in  the  stable  region. 
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TABLE  6-3 


VARIATION  OF  "A"  AND  "B"  VERSUS  "T"  FOR  CONSTANT  OPERATION 

ON  (un,0  =  (1.63. 8) 


T 
sec 

A 
units 

B 

units 

K 

P 

0 

0.0 

7340 

14680 

95 

190 

0.5 

2853 

2962 

37 

38.4 

0.1 

5722 

10568 

74.1 

136.6 

0.15 

5574 

9453 

72.2 

122.5 

0.2 

5280 

8279. 

4 

68.4 

107.2 

0.25 

4893 

7130 

63.3 

92.3 

0.3 

4463.6 

6064. 

1 

57.8 

78.6 

0.35 

4026.4 

5110. 

3 

52.2 

67.5 
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E.    TIME  DOMAIN  SIMULATION 

■  1.    It  would  be  convenient  to  present  at  this  stage  an  appropriate 
transformation  of  the  system  from  its  physical  block-diagram  to  another 
mathematically  equivalent  diagram  which  is  appropriate  for  simulation. 
Two  such  mathematical  block  diagrams  are  shown  in  figures  6-2a  and  6-2b. 
The  conventions  used  in  figure  6-2a  are  shown  in  table  6-4. 

2.  From  this  time-domain  simulation  the  time  responses  of  figures 
6-3a,  to  6-3f  were  obtained,  with  the  results  as  shown  tabulated  in  table 
6-5. 

3.  The  steady  state  output  one  gets  from  equation  (VI-1)  is  given 
by  the  formula 

1  *      36,7 


ss       .48313  -  147.59 

F.    CONCLUSIONS 

1.  The  system  proved  to  be  stable  for  T  -■  .1,  .15,  .2,  .25,  .35. 

2.  The  results  of  section  IV  and  the  time  simulation  coincide  as 
far  as  the  T '  =  . 1,  T=  .2  cases  are  concerned. 

3.  For  the  T  =  .3  stability  analysis,  it  seems  that  for  both  the 
section's  IV  method  and  the  time  domain  simulation,  one  would  expect  to 
get  stability  for  T  =  . 3  by  trying  bigger  frequency  and  final  time  ranges 
respectively.   But  since  it  was  decided  to  use  T  =  .15,  this  was  not 
necessary. 

4.  The  system  is  safe  stability-wise  to  use  for  the  selected 

parameters  of: 

A  =  5574      K  =  72.2 
P 

B  =  9453      K,  =  122.5 

0 
with  a  time-delay  constant  of  T  =  .15. 
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The  time-  domain  characteristics  for  a  step  response  are  for 

this  case: 

"Max.  Output"      _  - 
a*   Steady  State  "Output" 

P.   Time  for  "Max  Output"    =  9.19  sec 

Y.   Rise  Time  (107o-90%)     =  4.86  sec 

6.   Time  to  reach  (1%)  error=  10.4   sec 

e.   Error  at  "t=20  sec"     =  2.8   ppt 
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TABLE  6-4 


1. 


Delay 

-Ts 
e    = 

Tc 

,(s) 

Numerator 

100 

Call  I 

= 

.305K 
D 

= 

L  = 

7 

then: 

N2  =  L 

X 

Td(s 

) 

77.2,  (.444)  =  LL,  (.488)  =  L2 


N  =  L  x  L  x  T,(s) 

N  =  L.  x  L  x  T,(s) 
o    2        d 


3.    Denominator 


Dr    =  1.0 

6 

D5  =  98.469 

D,  -  4941.1680 
4 

D.  -  8978.082  +  A  x  T(s) 
j  d 

D_  =  5086.9961  -i-  [(.444)  x  A  +  b]  x  T,(s) 
2.  d 

Dn  =  3162.6011  +  [(.488)  x  A  +  (.444)  x  b]  x  T,(s) 
-1  d 

D  =  -147.59  +  [(.488)  x  B]  x  T,(s) 
o  d 

DD3  =  8978.082 
DD2  -  5086.9961 
DD1  =  3162.6011 
DD0  -  -147.59 
4.   With  the  conventions  of  table  6-2  at  hand,  the  system's  transfer 

function  becomes: 

N„(s)s   +  N.  (s)  s  +  N  (s) 
c  .    v     2   1   o   

R  6         5         4         3         2 

D6(s)s   +  D5(s)s   +  D4(s)s   +  D3(s)s   +  D2(s)s   +  D1(s)s  +  D  (s.) 
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or 


0 

R 


2   4  15  o      6 

s  s  s 


5s  42  33  24  15  o      o 

s  s  s  s  s 


From   this    form   the   signal- flow  diagram  of    figure   6-2a,    follows   easily, 
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APPENDIX  A 

ADDITIONAL  INFORMATION  CONCERNING  THE  ROOT-LOCUS  ANALYSIS 
OF  A  SYSTEM  WITH   TIME- DELAY 

This  appendix  contains  more  specific  information  concerning  the 
Pade   approximations  as  well  as  data  for  the  characteristic  equation  of 
the  system: 


C  =  N'  (s) 
R  =:  D '  (s  ) 


(see  Table  II) 


resulting  from  the  block  diagram   of  figure  2-2,  but  with  the  above 
mentioned  approximation  inserted. 

The  Pade  Approximation 

If  one  represents  by  DE(s),  the  transfer  function  of  a  time-delay 
in  the  "s-domain"  then  for  T  being  the  time-delay  constant: 


DE(s)  =  e~  !  =  1  -  Ts  + 


1MI  .  ll£L 

21      3! 


+ 


As  already  mentioned  in  section  II,  the  objective  of  the  Pade  ap- 
proximations is  to  define  two  polynomials  N  ,  (Ts)  and  D    (Ts)  of  order 

a ,  l>  a ,  d 

"a"  and  "b",  respectively,  which  have  a  quotient 


a,b 


rN    (Ts) 
a ,  b 

.Da  ,  (Ts)  _ 
a ,  b 


such  that  a  certain  number  of  terms  of  both  the  related  series  expansion 
of  DE(s)  and  this  quotient  are  identical. 
Letting: 

N  ,(Ts)    X      +   X,(Ts)  +  MTs)2+  ...  +  *■  (Ts)3 
F    (Ts)  L    a?b =  _o L_ ? a 

a'b   '      Da,b(Ts)   ^Q  +n1(Ts)  +^2(Ts)2+  ...  +Hb(Ts)b 
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one   observes   that   F        (Ts),    is   a    function  of  the  order   a,b   of  N        (Ts) 

a  j  d  a  j  b 

and  D  ,  (Ts).   So  after  selecting  a,b,  one  must  choose  the  appropriate 

ci  j  D 

values  for  the  coefficients  X     t  A.   and  u   t  Lk,  in  order  to  define  the 

o    a      o    b 

Pade  approximation. 

The  polynomials  N  K(Ts)  and  D    (Ts)  are  specified  by: 
a  9  D  a  y  o 


n=a  m=n~l 

n 
^)        ^  a-m 

N 


(Ts)    =    1  +V(-l5-^-     TT 


a  ,b  L,  n!  a   +  b-m 

n=l  m=0 

n=b  m=n-l 

D     K(Ts)   =   1  +yi~^       TT 


<,JLS;     =      1     +     /  j  II  " 

a,b  L,     n.  a  +  b-m 

n=l  m=0 


In  Table  Al ,  a  list  of  the  first  six  Pade  approximations  is  tabulated. 
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TABLE   Al 


LIST  OF   THE   FIRST   SIX  NORMAL    (a=b)    AND   REDUCED    (a=b-l)    APPROXIMATIONS 


Call   Ts   =   X: 


*i    i  (X)    -  T-T-f  Fn    ,  (X) 


1,1  x/         2   +  X  0,1  v    '    "    1   +   X 

6X  --  X2  „         „rX  6   -    2X 


F,    ?(X)    =  =-  F         (X)   = 

12  +  6X  +  X  l'  6  +  4X  +  X 

v        M        120   -    60X  +   12X2   -    X3  _        ,_        60   -    24X  +  3X2 

120  +  60X  +   12X     +  X  '  60  +  36X  +  9X     +  X 

™        1680   -   840X  +   180X2   -    20X3  +  X4 
F         (X)   =  "34 

1680  +  840X  +   180X     +  20X     +  X 

,„*        840   -    360X  +   6 OX2   -    4X3 

3   4W   "  "  "  2  3  4 

J'H  840  +  480X  +   120X     +   16X     +  X 

,-rt        30240   -    15120X  +  3360X2   -   420X3  +  30X4   -    X5 
F         (X)   =        ■  ^  3  4  5" 

30240  +   15120X  +  3360X     +  420X     +   30X     +  X 

.   .    _   15120   -    6720X  +   1260X2    -    120X3  +   5X4 
4'5  15120  +  8400X  +  2100X2  +  300X3  +  25X4  +  X5 

9  o  /,  c  €L 

665280  -  332640X  +  75600X  -  10080X  +  840X  -  42X  +  X 
665280  +  332640X  +  75600X^  +  10080X  +  840X  +  42X  +  X 

.   .        332640  -  151200X  +  30240X2  -  3360X3  +  210X4  -  6X5 

332640  +  181440X  +  45360X  +  6720X  +  630X  +  36X  +  X 


Note:   For  more  detail  on  the  Pade  approximation,  see  reference  [3J 
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TABLE  A2 

DATA  RELEVANT  TO  THE  CHARACTERISTIC  EQUATION  OF  THE  SYSTEM  SHOWN  IN  FIGURE 

2-2  OF  SECTION  II 

Characteristic    Equation: 


(0)s°  +    (1680)s1  + 


840T 

2  . 
s  + 

r      21 

180T 

s3  + 

3* 
20T 

4 

s   -t 

r  t*i 

s5  + 

+ 
1680 

+ 
840T 

+  2 
180T 
•    < 

+3 
20T  . 

T4 

*         * 

6 

+ 


k[(1680)s°   -    (840T)s  +    (180T2)s2   -    (20T3)s3  +    (T4)s     J      =      0 


Calculations  For  Repetitive  Terms 
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APPENDIX  B 


PARAMETER  PLANE  (CONSTANT  "x  "  and  "£")  EQUATIONS  FOR  A  SYSTEM  WITH  A 

TIME  DELAY 


The  characteristic  equation  of  a  system  is  given  to  be: 
n        i 

Z,VsK  +£a£  e'Ts  sK  =  0  (B-l) 

k=0      k=0 

where    a  *'=  Ab^  +  Be  (B-2) 

K.      K      K 

These  are  the  equations  (III-2)  and  (III-3)  which  characterize  the 

system  of  figure  3-1  encountered  in  chapter  3. 

The  objective  is  to  derive  a  set  of  two  equations  of  the  form: 

A'B.  (x  ,C)  +  B'C.  (;;  ,C)  +  D,  (•  ,D  =  0 
1   n  s        1   n '  ^      In' 

A*B2(0)n,C)  +  B.C2(xn,C)  +  n2(Xn,C)  -  0  (B-3) 

r  "i 

where:   i  B.  ,  C,  D.  should  hs   explicit  functions  of  X  ,  and  F 

L  l    l    lj .  n  „  n      s 

i=l,2 

and  A,B  the  unknowns. 

Using  the  fact  that 

i9 
s  =  x  (cos  8  +  i  sin  9)  =  x  eJ 
n  n 


it    follows    that: 
i 


,.      -Ts      K       V  ..      "TJU    (cos  8    +   j    sin  8),         i8 . 

LaKc      *s    ^LaKe     n  (V    > 

k=0  k=0 

I 

V  .,      -Tx   cos  8    -    jTx   sin  8        K     jK9 
ya      e        n  n  UJ        eJ 

LK  n 

k=0 

JL 

\    .,        K     -Txncos  8      j[K0    -    Txnsin  8] 
/■a!  x        e        n  eJ  n 

LK     n 

k=0 


K 
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»        a-  J^KG  -  TJU  sinG]      ,     _.     -Tincos  9       _  _ 

Now  expanding  eJ        n        and  getting  e    "      out  of  the 

summation  operator,  one  can  get: 

I 

£ajj  e"Ts  sK  =  (B-4) 


k=0 

i 


/a"  uu  Kl  cos(KB    -   TtDnsin  9)   +  jsin(K9    -   TUj  sin  9)1 


-T^ncos  i 

k=0 
Similarly   then: 

n  n 

)6V   sK  =  )d  ^  K(cosK9    +  jsin  K9 )  (B-5) 

Z_j  K      L-j   K  n 

k=0      k=0 
By  substituting  equations  (B-4),  (B-5)  into  equation  (B-l)  and  by  separa- 
ting the  equation  into  two  parts,  one  for  reals  and  one  for  imaginaries , 
it  follows: 

n  I 

)d,;x   K   cos   K9    +  G"TXncos        )a"    r   K  cos(K9    -    TV    sin  9)    =   0 
LKn  L  K     n  v  n 

k=0  0 

n  i 

)d  w  K   sin  KB    +  e~TinCOS   (     )  a"   UJ   K  sin(K9    -    Ttt   sin  9)   =   0 
LKn  ZL  K      n  n 

k=0  .  0 

and  by  substituting 

a"  =  Ab  +  Be  ,  (B-2),  one  arrives  to  the  following  two  equations 
K     K     K 


A  |yTxncos  9  ^K   cos  (K0  _  TjL^s.n  Q)j  + 


k=0 

i 


B  re"TiJncos  9  Vc  U)  K  cos    (KB    -    Ti^sin  9)1      + 


k=0 
n 


)d     x  K  cos   kg]   =   0  (B-6a) 


k=0 
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and 


I 

.   f    -Txn    cos  9   7",         K      .     ._         M        ,      ovl 
Ale        n  ybfJU        sin(K9    -    Tju   sin  9)      + 

L  L.  K  n  n  J 

k=0 

„  f    -TUUn   cos  9  V  K      .     ._         m        .      AN1    , 

Be        n  ,c^u        sin(K9    -    Tii   sin  9)      + 

k=0 

n 

}d  iju  K  sin  KD  j    =   0  (B-6b) 

L  Li  K  n  J 

k=0 

One  should  notice  that  equations  (B-6)  are  more  or  less  of  the 

form  of  equations  (B-3),  which  are  the  target,  since  cos9 ,  sin9 ,  and  9 

can  be  expressed  as  functions  of  £,  so  theoretically,  the  problem  is 

solved.   But  since  in  general  a  digital  computer  solution  is  desirable, 

it  is  noticed  that  a  summation  in  K  containing 

sin(K9  -  Ti'  sin  9)  or  cos  (K9  -  TJU  sin  9) 
n  n 

will  not  be  too  easy  to  evaluate. 

So  this  is  the  point  that  one  must  refer  to  reference  [l]   and 
appreciate  the  fact  that  he  found  some  existing  recursive  relations  in 
his  parameter  plane  analysis  without  a  time-delay.   From  this  fact,  it  is 
of  value  to  think  that  perhaps  the  same  thing  applies  for  the  case  under 
consideration 

For  this  reason  one  can  try: 

1.    First  Set  of  Substitutions. 

Call: 

Y  =  sin(K9  -  Tl  sin  9)  (B-7a) 

K  n 

X  =   cos(K9  -  TJJ  sin  9)  (B-7b) 

2h 

From  equation  (I-7b),  by  substituting   cos  9  =  -£ ,  sin  9  =(l-£  ),  and 

by  appropriate  trigonometric  expansion  it  follows  that: 
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X-l   ■   "(v''1-?2   Sin  b»nM2}   +  C   co.  [T«.n./l-C2  ]) 


X0     =c 


cs  [T,n  yi-c2] 


Xj   =  V'l-C2  sin  [t^  ./l-{2  j  -  {  cos  [xnn  Jl-Z2   j       (11-8) 

Similar  relations  can  be  derived  for  Y   ,  Y  ,  Y  but  this  is  not 
necessary. 

It  can  be  shown  by  induction  applied  to  equations  (B-7),  (see  Appen- 
dix C) ,  that: 

\+l  "  "2«\  "  *K-1  (B"9a) 

\+i  -  -2r\  -  Vi  (B"9b) 

\    -xatc^  +  vi"1  (B"9c> 

2 .         Second  Set  of  Substitutions. 

Call: 

U  -  sin  K9  (B-lOa) 

K 

T  =  cos  K9  (B-lOb) 

K 

It  follows  that: 

T-l  =  < 

To    =1 

Tj  =  -C  (B-ll) 

and  again  by  induction  from  equations  (I- 10)  it  can  be  proved  (see  Appen- 
dix C),  that  : 

TK+1  =  -2^TK  "  TK-1  (B"12a) 

UK+1  "  "2«UK  "  Vl  (B-12b) 
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Now  after  noting  that 

and  by  using  equations  (B-7)  and  (B-10)  into  equations  (B-6),  one  can 
get  equations  of  the  form  (B-3)  with: 

Bl  "  W  C'  V 

ci "  W  c  V 

Dl  "  W  «»  V 

B2  ■  W  c  V 

C2  ■  C2<V  C.  V 


D2  "  W  ^>  V 


(B-13) 

The  final  step  in  this  derivation  is  to  use  equations  (B-9c)  and  (B-12c) 

which  provide  Y  and  U  ,  in  terms  (X.,  X   ,)  and  (T  ,  T   .)  respectively 
K       k  K    K- 1         K    K-  i 

So  finally  from  (B-6) , (B-13) ,  (B-9c),  (B-12c)  it  follows: 

A  i  eT'^    >b  UJ  KX    (C>   «J    )]    +  B      eTJ^  )c.UJ   KX    (£ ,    U>)     j   + 

L  Z_i  K  n     K  n  J  L  Z..  k   n      K.  s         n     j 


K=0 


K=0 


[  ZV»\«>  ]   =  ° 


(B-14a) 


k=0 


^  ZVnVVCV  +  *k-l«-V)]  + 


K=0 

I 


B  [eT^?  Xv^e-V-V  +  *k-l«>.>)]  + 


K=0 


[  iVnK  (c"x<c>  +  Vi«>)]  "  ° 


(B-14b) 


K=0 
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where: 

TK+1   =   -2^TK  "    TK-1  (B"12a) 

hn  "  -2Z\  '  \-i  (B"9a) 

Equation    (B-8)    and    (B-ll),    giving   X      ,    X    ,    X     and    T      ,    T    ,    T     are 

also   necessary   if   a   computer   solution  of   equations    (B-14)    is    implemented. 

The   definitions    of   [B.,    C,    D.].    ,    0    follow   by   comparison  of   equa- 

1        1        1    i=l,2 

tions    (B-3)    and    (B-14): 

I 

tju  r  ci      k. 


\-^L\\X<Z'V 


K=0 

i 


Cl  "  ^  IVnV'V 


K=0 


n 

K_ 


D     =  ;d JU  "T   (£) 

1        uj  K   n      K 

K=0 


B2-e^Zv^«'V+Vi<t-  V] 


K=0 

K=0 
n 

D2=Id^nK[«TK«)+  Vl«)]  (B'l5) 

K=0 

So    if   B   C      -    B   C     r   0   one   can   get    from    (A-3): 

A-    (CXD2    -    C^J/CB^    -    B2Cl) 

B  =  (B2D1  -  D2B1)/(B1C2  -  B  C^  (B-16) 
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APPENDIX  C 
PROOF  OF  THE  VALIDITY  OF  THE  RECURSION  RELATIONS  USED  IN  APPENDIX  A 

1.  To  prove  that  the  relations  (1-9)  are  true  one  must  first 
prove  their  validity  for  K  =  1,  and  then  given  the  truth  of  these  rela- 
tions for  K  =  V  one  must  prove  their  truth  for  K  =  V+l . 

2.  Validity  of  Relations  (B-9)  for  K  =  1. 

a.    In  that  case  it  is  required  to  prove  that: 

(1)  X2  =  -2CXj_  -  XQ  (C-l) 

(2)  Y2  =  -2CYX  -  YQ  (C-2) 

(3)  Yj_  =  (1-C2)"%[CX1  +  XQ]  (C-3) 

given,  of  course  that, 


and 


Y„  =  sin  (KB  -  Tx  sin  9) 
K  n 

Y  =  cos  (KB  -  Td  sin  9) 

b.    By  direct  substitution  of  X  ,  X  ,  X  and  appropriate 
trigonometric  expansions  in  equation  (C-l)  one  gets: 

cos  29  cos(T.f  sinB)  +  sin  28  sin(Tx  sin  B)  = 
n  n 

(-2£cos8  -  1)  cos  (Til!  sinB)  -  2£sin9  sin(Tx  sin  6) 

and  taking  into  account  that: 

C  =  -cosB 
it  is  obvious  that: 

-2£  cos  0  -1  =  2  cos  B  -  1  =  cos  28 
and  -2£  sin  8=2  sin  8  cos  8  =  sin  28 
Therefore,  equation  (C-l)  is  true. 
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c.  In  a   similar  way  one  can  get   that   equation    (C-2)    is    true. 

d.  Substituting    for   Y.,    X    ,    X      in   equation    (C-3),    one   gets    the 
following: 

sin  6    cos  (TJU   sin9)    -   cos  9    sin(TJJ   sin9)    = 
n  v      n  ' 

(1-C    )    2[Cc°s9  •cos(TX'   sin9)    +  £sin9  «sin(TJU   sin9  )+cos  (-Tx   sin9)] 

ny^         vn/vn 

This  relation  is  true  since  £  =  -cos  9  as  mentioned  before.   So  equation 

(C-3)  is  valid. 

3.  Now  in  the  case  that  K=V,  one  must  prove  that: 

(1)  XV+2  =  -2^Vl  "  XV  (C"7) 

(2)  Yv+2  =  -2CYv+1  -  Yy  (C-8) 

(3)  Yv+1  =  (1-C2)"%[£X1  +  XQ]  (C-9) 
provided  that 

(4)  Xv+1  =   -2CXv  -    XV_L  (C-4) 

(5)   Vi  =  -2^v  '  Vi  (c"5) 

(6)      Yv       -    (l-cV^CC^  +  Vl]  (C"6) 

4.  The  proof  for  the  three  sets  of  equations  in  pairs  of  two 
follows  the  same  pattern,  so  it  is  enough  to  actually  prove  that  the  re- 
lations (C-4)  and  (C-7)  are  true.   This  can  be  done  in  steps  as  follows: 

a.    By  trigonometric  expansion  of  the  relation 
X  =  cos  (K9  -  [Tjj  sin9]) 
for:  K  =  V-l,  V,  V  +  1,  V  +  2 

X^   =  (cosV9  cos9  +  sinV9  sin9)  cos[Tju  sin9]  + 

(sinV9  cos9  -  sin9  cosV9)  sin[Tx  sin9] 

n 

X  =  cosV9    cos[Tx   sin9]   +  sinV9    sin[Tx   sin9] 

V  n  n 

=  (cos(V+l)9    cos9    +   sin(V-+l)9    sin9)cos[Tx   sin9]   + 

n 

(sin(V+l)9    cos9    -    cos  (V+l)9    sin9)sin[Tx  sin9] 

n 
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X         =   cos(V+l)9    cos[TJJ   sin9]   +  sin(V+l)8    sin[TJJ   sin9] 
V+l  n  rt 

=    (cosVG    cos9    -    sinVe    sin9)cos[TJ0  sin9]   + 

n 

(sinV9    cos9   +  cosV9    sin9)sin[Tx   sin9] 

n 

*V+2  =    (cos(V+L)G    cos[Ti   sin9]   +   sin(V+l)9    sin[TU)   sin0])cos9    + 
(cos(V+l)9    sin[Tx   sinG]    -    sin(V+l)9    cos[Ti/   sin9])sin9 

b.  By  observing   that: 

X  =   X  cos9    -    sin9    sin[v9    -    Tju   sin9],    and 

I,   ,    =   X  cos9   +   sin9    sin[\6    -    Tju   sin9] 
v-1  V  n 

which   added    together,    they  give 

Vi  +  Vi  - 2  cose  V  or 
Vi =  "2^xv "  Vi 

X      9    =   X        cos9    -    sin9    sin[  (V+l)9    -    TiJ   sin9]      and 
X  =   X       cos9   +   sin9    sin[  (V+l)9    -    Td   sin9] 

which   added    together  give 

V2  +  Xv=  2  cose  Vi'  or 

\+2  =   "2CVi   "    *V 

c.  Now  since   Xr  -  =  -2^X      -    X    holds  for  K  =  1  as  well  as 
for  K  =  V  ,  and  V+l,  it  holds  for  every  other  value  of  K.   Q.E.D. 

5.    The  proof  that  equations  (1-12)  are  also  valid   follows  the 
same  procedure  discussed  from  paragraph  1  up  to  4  of  this  Appendix.   But 
in  order  to  prevent  unnecessary  duplication  one  can  show  only  some  portion 
of  the  last  part  of  this  proof  just  for  demonstration: 

a.    Since  Tr  =  cos  K9  and  U  =  sin  K9  by  trigonometric  ex- 


pansion  follows: 


(1)  T     ,  =  cosK9  cos9  +  sinK9  sin9 

T.   =  cosK9 
K 

T  =   cosK9    cos9    -    sinK9    sin9 

K+l 
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can  get 


(2)   U   .  =  sinKB  cos9  -  cosK9  sin9 
K-  1 

U,.   =  sinKB 
K 

U^  1  =  sinK9  cos9  +  cosKD  sin9 

b.    By  considering  the  expressions  for  T  ,  T   ,  and  U  ,  one 

K         K~  1  K 


Tv   i    =   T  cos9    +  Usin9,    or 
K-  IK  K 


=    (sin9)      (-cos9TK  +  TK_L),    or 


U 
K 

UK       =    (1"C2)'%(CTK  +  TK-1)  (C-10) 


c.        One   can  observe  also   that: 

T     _    =   cos9(   cosK9    +  -^|   sin  KB) 
K-l  cos9 

T„   ,    =   cos9  (  cosKB    -    S1"     sin  K9 ) 
K+l  cos9 

Therefore,    by   adding    these   two   equations,    the   recurrence   relation   needed, 

follows : 

TK+1   +  TK-1   "   2   C0S   K9    "   2TK'    °r 

TK+1   "   -2?TK  "    TK-1  (C"U) 

uk-i  =  coe8(uk-^cos,S) 
UK+1=.cose(UK  +  -|i|cosKB) 

and   again  by   addition   it    follows: 

\+l  "   "2«UK  "   Vl  (C"12) 
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APPENDIX  E 


SOLUTION  TO  A  QUARTIC  EQUATION 


1.    A  quartic  equation 


/.  Q  O 

X  +  aX  +  bX  +  cX  +  d  (E-l) 

has  the  resolvent  cubic  equation 

Y3  -  bY2  +  (ac-4d)Y  -  a2d  +  4bd  -  c2  =  0  (E-2) 

2.  By  using  the  attached  solution  to  a  cubic  equation,  one  can  specify 
arbitrarily  a  solution  y  =  y  ,  to  equation  (E-2)  and  then  define 

»  =  ■k  ■   b  +  yr 

3.  The  four  roots  to  (E-l),  will  then  be: 

_  •  a   R   D 

x  "  4  +  2  ±   2 

a   R  ,  E 
x  ~  "  4  "  2  "  2 

4.  The  parameters   E  and   D,    take   the   following  values: 

a.        R  ^   0 

— 2 3— 

/3a  2  4ab-8c-a 

D  =  ./— ; R  ~2b   +     — 

\    4  4R 


tt-2- 


/3a        „2      0^      4ab-8c-a' 


E  =  j=^-  -R    -2b- 


4R 
b.         R  =   0 


D  =  Pt~  -   2b  +  2/7~2    -   4d 
v    4  r 


rf 


E  =      ■%-  -   2b-    2/7"2   -   4d 
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ATTACHMENT  TO  APPENDIX  E 


SOLUTION  TO  A  CUBIC  EQUATION 


1,    Consider  the  Cubic: 

3     2 
Y  +pY  +  qY  +  r  =  0  (E-l) 

Let:  y  =  x  -  £.  (E-2) 

Calculate: 

2 

2  2   2px    p 
y  =  x  -  -f-     +  f- 

3  3    2      2    3 
y  =  x  -  x  p  +  xp  -  _p_ 

3    27 

and  substitute  into  (E-l)  ,  it  follows: 
2      3 
x  +  x(q  -  -2-0  +  -£-   "  "f9  +  r 

2 o    The  cubic  is  now  in  the  form: 

2  3 

3         ,  p       ,  ,    2p     pq 

x  +  ax  +  b  =  0,  where  a  =  q  -  -r-  ,  and  b  =  -rb-  —   -^  +  r. 

3.    Letting: 


—2 3—  —2 3— 

3/  b    /b    a         „  3  /   b    /b    a 


*     T   b  a.  J°i  B  M   3  /    b     / 

A°V"  2  +  vT  +27    >    M"V  "  2  "V' 


4  +  27 


The  three  solutions  for  y  come  out  to  be: 
yx  =  A  +  M  -  % 

A+M   A-M/- r 


A+M   A-M/-— 
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